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ABSTRACT 


Inviscid hypersonic flow of chemically relaxing air about pointed circular 
cones has been calculated using Dorodnitsyn's one-strip integral method. The air 
model used consists of C^, O, N, and NO, of which the relative amounts 
are controlled by six reaction equations which actually represent eighteen chem- 
ical reactions. Also considered are the limiting cases of nonequilibrium flow, 
namely, chemically frozen flow, and the flow with chemical equilibrium. Vibra- 
tional equilibrium has been assumed throughout. 

The principal contribution of this work lies in the method of solution. 
Specifically, it is demonstrated that the replacement of the approximate tangential 
momentum equation, and of all approximate species continuity equations by their 
exact forms along the body surface leads to a system of equations which, first, can 
be stably integrated, and second, yields results which compare quite well with 
those obtained by the method of characteristics. A detailed discussion of the 
polynomial approximations employed in the integral method, the introduction of a 
characteristic relaxation length defined in terms of the initial species mass fraction 
gradients, as well as many new results for hypersonic flow of a ; r with finite rate 
nonequilibrium dissociation about circular cones are further contributions of the 
work reported herein. 
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The discussion cf results includes two cases which show the very good 
agreement of the results from the integral method with those from the method of 
characteristics. Various other cases serve to illustrate the influence of free stream 
velocity (5400 2 Cm/sec] ^11186), cone semi -vertex angle (20° s 0 ^45 °), 
and altitude (31 km to 80 km) on the variation of all flow field variables. The 
analysis of a cone with constant semi -vertex angle and free stream velocity at 
three different altitudes is used to demonstrate the validity, and its range, of the 
nonequilibrium scaling law P M * L “constant. 

For one case, the asymptotic equilibrium state of the nonequilibrium flow is 
compared with the. results of a conical equilibrium flow calculation. According 
to theoretical considerations an agreement between the two cannot be expected, 
and this is confirmed by the comparison. Although the cone surface pressures 
agree very closely, it is found that the values for the asymptotic body surface 
temperature are higher than those for conical equilibrium. Both, the asymptotic 
surface density and the surface velocity stay below their values for conical 
equilibrium . 
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CHAPTER I 


INTRODUCTION 

Over the past few years considerable effort has been devoted to the under- 
standing and the theoretical prediction of the nonequilibrium flow field conditions 
encountered in hypersonic flight. Present technology dictated the use of blunted 
body shapes for reentry vehicles in order to minimize the convective heating rate 
to the body surface. Quite naturally/ calculations of flow fields surrounding 
hypersonic blunt bodies were of prime importance. For a survey of calculation 
techniques for this class of problems the reader is referred to Hayes and Probstein 
(Ref. 1). In particular, the general method of integral relations, proposed by 
Dorodnitsyn (Ref. 2), proved to be quite successful in predicting shock wave shapes 
and other flow field variables for certain simple blunt bodies. Also, this method 
is readily adaptable to machine calculations. The investigations of Shi h, Baron, 
Krupp and Towle (Ref. 3), Hermann and Thoenes (Ref. 4), and most recently, 
Thompson (Ref. 5), are some examples of Inviscid real gas flow field calculations 
using the integral method. 

Comparatively less attention has been paid to the high temperature real gas 
flow past pointed bodies, which appear to be of growing interest and importance. 
Certainly the time will come when it is desirable to land on the Earth's surface 
some space probes returning from far distant planets. These probes will reenter the 
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atmosphere at a speed far above earth parabolic speed. Both, convective and 
radiative hea''ng will occur, but because of the high speed, radiation is expected 
to be the dominant source of heating . It was therefore pointed out by Allen (Ref .6) 
that the shock wave sweepback provided by pointed conical entry bodies could be 
used to reduce the gas luminosity grossly, and thus the total heating rate. There- 
fore, besides being of fundamental interest, the calculation of high temperature 
nonequilibrium flow past pointed cones may be of very practical importance in the 
near future . 

It is interesting to note that one of the first approximate theories for non- 
equilibrium cone flow was developed by Chapman (reported by Stephenson, Ref .7) 
in order to determine relaxation times of oxygen and nitrogen vibrations. Using 
essentially a stream tube method. Chapman derived equations which related the 
coordinates of the curved shock wave to an effective relaxation time of the flow 
in the shock layer. Having obtained shock wave coordinates from experiments, 
Stephenson (Ref. 7) evaluated relaxation times which agreed reasonably well with 
data obtained from shock tube experiments. 

If the flow field is entirely supersonic then the method of characteristics 
is available for its calculation, provided that some initial data are available. 
Calculations using the method of characteristics were first reported by Sedney, 

South and Gerber (Ref. 8) for the case of vibrationally relaxing flow of pure nit- 
rogen over a wedge. These calculations were later extended to include flows past 
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circular cones (Ref. 9). The authors also reported the necessity of introducing 
modifications in the method of characteristics in order to increase the accuracy of 
the results. This work was followed by Capiaux and Washington (Ref. 10) who 
treated nonequilibrium flow of an "ideal dissociating gas ' (so-called Lighthill gas) 
past a wedge, also using the method of characteristics. 

In Ref. 11, Wood, Springfield and Pallone investigated the flow over axi- 
symmetric bodies by using the method of characteristics. They employed a multi- 
component gas model and included a simple model for vibration-dissociation 
coupling . Results were given for flows over blunted cones and ogival bodies. 

The first application of Dorodni tsyn ‘s integral method to the real gas flow 
over pointed bodies was given by South (Ref. 12). He considered vibrational ly 
relaxing flow of a pure diatomic gas past a circular cone, using the one strip 
approximation, and constant relaxation time. He also reported that a system of 
equations containing the tangential momentum equation and the vibrational rate 
equation in their exact forms could not be integrated successfully. He found it 
necessary to use linear approximations in all the conservation equations. The 
pressure distribution obtained from this set, together with a set of corrected 
equations then served to obtain a corrected solution. The application of the two 
and three strip integral method to the same problem is discussed in Ref. 13. 

South and Newman (Ref. 14), and Newman (Ref. 15) discussed a modifi- 
cation to the one strip integral method as applied to wedge flow. Instead of using 
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first order polynomials to approximate certain integrands, the integrals themselves 
were replaced by weighted averages of the integrands* According to the authors, 
this procedure resulted in a better description of the flow within the shock layer, 
but again, it was found necessary to integrate a set of corrected equations after 
obtaining an approximate solution. 

Thoenes (Ref. 16) used the one strip integral method and a simplified three 
component air model, applicable in the oxygen dissociation regime, to calculate 
chemically relaxing flow past cones. In particular, these calculations successfully 
employed the exact tangential momentum and species continuity equation, a method 
which the author labelled as the "semi -exact procedure". Results from the two 
strip integral method for dissociating flow were discussed in Ref. 17 . 

There are other methods which have been used for the calculation of non- 
equilibrium flow fields around pointed bodies- Lee (Ref. 18) evaluated vibrationally 
relaxing flow over wedges by a perturbation method, while DeJarnette (Ref. 19) 
presented an "artificial viscosity" method to evaluate wedge flows with simple 
relaxation processes, including various schemes of vibration-dissociation coupling. 
All these calculations with simple rate processes are of basic interest and valuable 
for determining suitable numerical techniques. However, with the exception of 
Ref. 11, they cannot describe adequately a flow field involving several simul- 
taneous relaxation processes as they occur, in air for example. Presumably, this 
prompted Spurk, Gerber and Sedney (Ref. 20) to calculate hypersonic flow of air 
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past wedges and cones by the method of characteristics, considering a five compo- 
nent air model involving numerous chemical reactions. Unfortunately, the appli- 
cation of the method of characteristics to nonequilibrium flow is complicated and 
computation times can be very long, even on modern computers. 

No calculations, comparable to the work of Spurk, Gerber and Sedney, but 
using procedures other than the method of characteristics, could be found in the 
literature at the present time. It is expected that the work herein will fill this 
gap by presenting a formulation of Dorodnitsyn's integral method for hypersonic 
flow of a multicomponent chemically reacting gas about pointed bodies. The 
principal advantage of the integral method is that it reduces the governing partial 
differential equations to ordinary differential equations for which highly efficient 
numerical integration techniques have been developed. 

A detailed description of the multicomponent gas model for high temperature 
air, and its thermodynamic properties, is given in Chapter II. In Chapter III, the 
basic partial differential equations are presented. They are cast in a suitable 
coordinate system, and the discussion of applicable boundary conditions concludes 
the formulation of the problem. The solution is discussed in Chapter IV* After 
reviewing the method of integral relations in general, its application to the 
problem at hand is described in detail • A final section on the numerical techniques 
used leads to Chapter V, which is aevoted to the discussion of the results. 



CHAPTER II 


GAS MODEL AND THERMODYNAMICS 

2 . 1 Gos Model 

To determine the values of the flow variables within the shock layer ooout 
a body moving in the atmosphere at hypersonic speeds/ it is necessary to account 
for the effects of chemical reactions that take place within the layer. Such 
reactions proceed at a finite rate/ but for flow about infinite cones/ the gas will 
eventually reach a state of thermal and chemical equilibrium. Figure 1 presents 
lines of constant compressibility factor and temperature at the surface of a cone 
for equilibrium flow/ using data given in Ref. 21 . Although this is a purely hypo- 
thetical case, the figure shows approximately the conditions which have to be met 
by the air model to be chosen. For example/ if a conical vehicle with a semi- 
vertex angle of thirty degrees enters the Earth's atmosphere at an altitude of 80 km 
at earth parabolic speed it is to be expected that/ once equilibrium is reached on 
the streamline that wets the body, the surface temperature is around five thousand 
degrees Kelvin, and that a sizable fraction of the nitrogen molecules must have 
dissociated (Z^l-4). This example shows that, in order to calculate such flow 
conditions, the gas model representin'! high temperature air must be assumed to 
consist at least of a mixture of C^, NO, O, and N. It is fjrthermore 
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assumed that all species are in translational, rotational and vibrational equilibrium 
in every flow region, which implies a single temperature for oil the degrees of 
freedom. Assuming weakly interacting particles (I .e. particles which exchange 
energy only by collisions, and the collision times are small compared with the time 
between collisions) partition functions may be used to calculate all thermodynamic 
properties. Thus ionization, vibrational nonequilibrium, as well as vibration- 
dissociation coupling are neglected. The latter two assumptions will be critically 
examined at the end of this chapter. 

The composition of air along with other physical constants as used in the 
present investigation are summarized in Table I. 

2.2 Partition Functions 

Under the assumptions stated in the previous section, all of the thermo- 
dynamic properties of a gas may be derived from its partition function. Details 
on the partition function may be found in any text on statistical thermodynamics 
(Ref. 22, 23, 24). In this section only those relations will be summarized which 
are needed for subsequent calculations* 

The partition function may be defined as 

€ 

o = ]£ 9j exp (- Tt- ) (2J) 

j 

where e. is the common energy of several states, and g. is the degeneracy, that 
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is, the nurrber of states of the particle which have this common energy level . In 
the present investigation, the energy may be due to the translational, rotational 
or vibrational motion of the particles and, in a limited way, due to the motion of 
the electrons within the particle. Assuming that no coupling exists between the 
different modes of excitation, the partition function may be written as the product 

Q=Q t Q r Q v Q e (2-2) 

where the factors on the right hand side of Eq. (2.2) are the partition functions 
associated with the translational, rotational, vibrational and electronic energy 
levels of the particle. For diatomic molecules these factors are: 


o f = v ( 2 nm 2 k t y /2 

(2.3) 

_ r 8 n 2 I k T _ T 

Q ^ - 

H h 2 K 0 r 

(2.4) 

Q V = £ 1 - exp ( - 0 V / T ) J 

(2.5) 

n 

Q 6 = £ 9: exp (- 0? /T ) 

j ”0 

(2.6) 


Here 0 1 ", 0 V , and 0? are characteristic temperatures for rotation, vibration, and 
electronic excitation, respectively. Monatomic particles have no modes of rota- 
tional or vibrational excitation, therefore, the respective partition functions take 
the value unity. Table II presents the atomic and molecular constants which were 
used in the present calculations. 
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2.3 Thermodynamic Properties 

2.3.1 Thermal Equation of Stat e 

According to statistical thermodynamics the partial pressure which a system 
of weakly interacting particles of species i exerts on its surroundings is given in 
terms of the partition function as 


p. = n, H -rT 7 “ (In O) (2*7) 

i i o V 

where n. denotes the number of particles of species i contained in the volume 
i 

V, and k is Boltzmann's constant. Since, as seen from Eq • (2.3) through (2-6), 
only the translational partition function depends on the volume V, one obtains 


W-('"Q) “ TT^< lnQt > 


I 

V 


(2.8) 


Hence, the pressure becomes 


n. k T 

t 


(2.9) 


Introducing the molecular weight M., and A vogadro's number N*, Eq* (2.9) can 
be rewritten as 

which is at once recognized as the equation of state of a thermally perfect gas, viz 



T 


(2.H) 



Since it is assumed that each species of the gas mixture can be treated as a 
thermally perfect gas, the pressure of the mixture is, according to Dalton's lav/, 
the sum of its partial pressures, that is 



Defining the mass fraction of the ith species as 

P. 

- _ i 


the equation of state can be expressed as 


( 2 . 12 ) 


(2-13) 


p * p R* T E Y | 
i 

where the Y. are the molar concentrations defined as 

i 



(2.14) 


(2.15) 


It is often convenient to express the pressure of the dissociated gas mixture in terms 
of the gas constant of the undissociated gas. Eq. (2.14) can then be wiitten as 


p = p RZT 


(2.16) 


w h e re R » R* /M is the specific gas constant of the undissociated gas, and Z is 
the compressibility factor defined by 
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A gas mixture consisting of i distinct species in the dissociated state, and 
of i species in the original undissociated state, requires i Q mass balance 
equations and (i - i^) relations between the species mass fractions. The former 
are, first, that the sum of all mass fractions is unity, 

£C. = 1 (2.18) 

I 1 

and, second, that the mass fraction of one element is fixed. For example, for air 

C C + C Q2 + (“^ 0 ) C NO = < C 02>o < 2 ' ,9 > 

or, alternatively 

C N + C N2 + ("M^O ) Cn ° ” ^ C N2 ) o * 2 * 20 * 

which is a direct consequence of relations (2.18) and (2.19). 


2.3.2 Internal Energy and Enthalpy 

Again calling on statistical thermodynamics, the internal energy of a system 
of n. weakly interacting particles is given in terms of the partition function as 


E. * n. k T 2 ' (In Q) (2-21) 

Introducing the molecular weight and Avogadro's number, the internal energy 
becomes 


E. 

i 


n 


M. 


N* k 

i 



b 

dT 


N* 


(In Q) 


( 2 . 22 ) 
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Recognizing that the term (n. M. / N*) has the dimension of mass, the specific 


i i 


interna! energy, that is internal energy per unit mass, can be expressed as 


= R* .23 

e . i j T i t (In Q) 


i 


M. 


i 


3 T 


(2.23) 


Using Eq. (2.2) through (2-6), and (2.23), the specific internal energy fora pure 


monatomic gas is 


e . 


R* ( 

M. \ 


T + 


S g . 9. 6 exp ( - 9? / T) 
J J J J 1 

£ g exp ( - 9.® / T ) 

j J J 


I 


(2.24) 


and, for a pure diatomic gas. 


e. 


- J$1 

M 


?{*- 


? g.0? exp (-0.® / T) 


T + 


+ -L-J-1 


exp (9. V /T) “1 


£g exp (~9. e / T) 
J J J 


! 


(2.25) 

The last terms in the above equations, representing the contribution of electronic 
excitation to the internal energy, become significant at extreme temperatures only, 
and can be neglected for the range of temperatures considered in the present 
investigation. 

In the case of a dissociating gas mixture the dissociation energy for each 
molecular species must be included in the internal energy. This energy term may 


be expressed as 


d _ 


e. 

J 


_R* 

M 


- r (C.) -C.l 9* 

j L j'o J J J 

J 


(2.26) 
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where, in the present case, j “■ N 2 , NO, and the subscript o designates 

the undissociated state. 

The total specific internal energy of a mixture of monatomic and diatomic 
gases can now be written as 



The specific enthalpy of the mixture is defined as 


(2.27) 


h = E + (2-28) 

P 

and substitution of the pertinent relations into Eq. (2.28) yields, for the five com- 
ponent gas mixture considered herein, the following expression for the enthalpy: 


h 

R *|t T(Y 02 + Y N2 

+ y no> 

+ 4 

- T (Y o +y n ) 



Y fl V 

, 02 02 

v 

, Y N2 

0 N2 

+ 

< 

z 

O 

CD 

z < 
O 



9 02* 
e - 1 

© 

CD 

2 < 
KJ 

— i 

1 

9 no /t , 

e - i 


+ [^02^0 

- Y 1 9* + 

t 02 J 02 

[^N2^o 

" Y N2l 

* 

0 N2 

+ [^NO^o " Y No] 








(2.29) 


Utilizing the mass balcnce equations, Eq. (2- 19) and (2.20),in order to eliminate 
the molar concentrations (Y Q2 ) o and (Y N2 ) q , and also considering that in the 
adopted air model (Y N q) q is zero, a more convenient expression for the enthalpy 
is easily derived from Eq. (2.29): 
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h = 


{ 


R* { Y 


02 1 2 


T + 


02 


9 V /T 

W A^ / 9 


+ Y 


N2 


fl v 

f T + — ^ 


02 


-1 


0 V / J 

N2 / . 

e -1 


+ Y 


NO y 


/ 7 


T + 


NO 


02 


N2 


6 nV . 

e -1 


- 0 


NO 


+ V. 


i_ T + -®°i 


+ y 


N 


i_ T +!n1 

2 2 


(2.30) 


If is recognized that fhe specific enthalpy now has the more general form 

h = E C. h. (I) 

: i i 


(2.31) 


For the calculations in the following chapters it is convenient to define 
a frozen specific heat as the partial derivative of h with respect to T, that is 


dh 

dT 


which, according to Eq. (2.31) may be expressed as 


(2.32) 


d h. 


E w ri. 

C. J = E C. c 

i o T ; i p. 


(2.33) 


p i o i j i P 

Hence, like the specific enthalpy, the frozen specific heat of the gas mixture is 
simply the weighted sum of the individual specific heats. An explicit expression 


for c^, as derived from Eq. (2.30), is given in Appendix A. 
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2.4 Rate Processes 

2.4> 1 Finite Rote Nonequilibrium 

In general/ a large number of collisions among the particles is required to 
equilibrate the molecular vibrations/ dissociation and higher modes of excitation 
with the local translational temperature. This means that, once being in a state 
of nonequilibrium/ a finite amount of time is needed for the gas properties to 
achieve thermodynamic equilibrium. The departure from equilibrium of a flowing 
gas is characterized by the magnitude of this relaxation time T relative to some 
translational time t needed by the particles to change their state- The limiting 
values of the ratio of t over t are a convenient means to characterize the 
equilibrium state (t / t - 0 ) and the frozen state ( t / t - 00 ). 

The purpose of this section is to present the equations which are necessary 
to account for finite rate nonequilibrium effects in the calculations. Essentially 
this is equivalent to finding an expression for the relaxation time T in terms of 
the variables of state. 

2.4.2 Species Continuity Equation 

Consider an arbitrary material volume containing particles of various species, 
each species forming the partial density p . Neglecting diffusion, the mass rate 
of change of the ith component in the mixture must equal the rate of production 
(which may be negative or positive) of the ith species, that is 


16 


jPt/p. dv = /w. dv 


(2.34) 


V V 

where W. is defined as the mass rate of formation of species i per unit volume 
due to chemical reactions. Using Reynolds' transport theorem, Eq. (2.34) maybe 


rewritten 


/[£ 


+ 7 • ( p. q ) 


-*■] 


dV - 0 


(2.35) 


Since the volume considered was arbitrary, this means that 

dp. ^ . 

~~ + ? • (p. q ) = W. (2-36) 

Summing overall species i, and considering that ^ p. - p, al so Ew. "0 

i 1 i 1 

since total mass must be conserved, one obtains the overall mass conservation 


equation 


dp 


at 


+ 7 • ( Pq ) = 0 


(2.37) 


Replacing p. by p C. in Eq. (2 .36), and expanding differentials yields 


c i [4r" + v ' (p ^] + p [ 


d C. 


i 


— + ^ 


7 C 


.] 


w. 


(2.38) 


Since the first term on the left hand side vanishes according to Eq. (2-37), the 


species continuity equation finally assumes the form 


dC. 

i 

at 


+ q-7C. “ (W. / p) 


(2.39) 


It is worth noting that the quantity ( p/W. ) has the dimension of time, and 
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therefore may be interpreted as a relaxation time for the species i . In the 
following section an explicit expression for W. will be derived. 


2.4.3 Rate Equations 

Following Vincent? and Kruger (Ref. 24), a chemical reaction is considered 


to have the following description: 


? V | A l +V M A M 


d,M 


E v.'A, + A,, 

T i i MM 


(2.40) 


Y,M 


where the v, are the stoichiometric coefficients of the reactants, those of the 

i « 

products. A is a so-called third body species which does not participate as a 
M 

reactant. k, tl and k are the respective reaction rate constants for disso- 
d,M r,M 

ciation and recombination. 

Defining the concentration of the species A. by X., that is 


c. 

X. = “TT““ o - P Y. 
i M. i 


(2.41) 


it is an experimentally observed fact that the rate of formation of any one of the 
products in a single step reaction can be expressed as 


v v 

i M 


d X. , 

~dT~) d = ( V i'V k d,M ^ X i X M 


(2.42) 


where the product on the right hand side extends over all the reactants. 
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Similarly, the rate of disappearance of any species A. can be expressed 


as 


/ d X. , v .' 

bf) r - ( V V r> fc r f M 1 ?^ x 


M 


M 


(2.43) 


where again the product is taken over all of the reactants- 

The net rate of production of the species A . is then the sum of its rote of 
formation and its rate of disappearance, vis. 


15.r d A) wli) 

dt \ d t / , \ d t / 


(2.44) 


or, with Eq. (2.42) and (2.43): 


dX. r 

-df = < v ;-vK 


M 


V. 


V 


M 


M V ^ X i'- k r,M V^V] 

(2.45) 


For thermodynamic equilibrium, the net rate of production of each and all 
species must vanish ("principle of detailed balancing"), hence from Eq. (2.45): 


v. 


k, xa n x. 1 - k .. jj x. 1 = o 

d,M j i r,M ; i 


(2.46) 


Eq. (2.46) can also be written in the more familiar form 


‘d,M _ 


n 


v 


X. 


‘r,W 


- K C (T) 


(2.47) 


IT x. 

i i 


where K (T) is called the concentration equilibrium constant, and expression 


(2.47) represents the law of mai,s action 
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If It Is assumed that, for given temperature and reactant concentrations, the 

rate at which a reaction proceeds in a given direction is the same whether the 

system is or is not in equilibrium, Eq. (2.47) may be used to replace either k^ ^ 

or k in Eq. (2.45). Hence let 
r,M n 


Using this expression, 
acting as third bodies. 


k , = k K (T ) 

d,M r,M c 


(2.48) 


and also accounting for the possibility of several species 
the resulting rate of production is obtained from expression 


(2.45) as 


d X. 


li- = (v.' - v. )J2 k I* X., M [ K Ilx.'-Xlx. '1(2 

t v i t ' r,M M L c : i \ i J 
M 


,49) 


Considering also that there may be m different chemical reactions of the type 
represented by Eq. (2.40), the total rate of production of the element A. is given 
by the sum over all reactions, that is 


d X. 

i 

IT" 


I -E 

tot m 


M 


v, 


X 


M 


M 


[ K TI 
L c j 


v. 


X 


“IT X. V i 


] 


(2.50) 

It is obvious from Eq. (2.47) that K (T) is, for the mth reaction, indepen- 

c 

dent of the third body species A . . . Naturally, for different reactions, the equili 
brium constants are different functions of temperature. Their evaluation is dis- 
cussed in section 2.5 of this chapter. 

For the air model used in this investigation the following set of chemical 
reactions will be considered (Ref. 15): 
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I 

k d,M 


I 

O 0 + M + 5. 1 ev £ 

2 k I 

r,M 

2 0+ M 

(2.51) 

II 

N 2 +M +9.8 ev ? 

2 N + M 

(2.52) 

III 

NO +M +6.5 ev ? 

O + N + M 

(2.53) 

IV 

NO +0 + 1.4 ev 

o 2 + n 

(2-54) 

V 

N 2 + O + 3.3 ev 2 

NO +N 

(2.55) 

VI 

°2 + N 2 + 1,9 eV - 

2 NO 

(2.56) 


Reactions I, II/ and III are the principal dissociation reactions, and, considering 
the various third body species (see Table III), actually represent fifteen different 
reactions. The remaining equations represent the so called shuffle reactions. 
Altogether, eighteen chemical reactions are considered between the five components 
of the gas mixture- 

The application of Eq. (2.50) to the above set of reactions yields the 
following expressions, already simplified by introducing the molar concentrations 

Y. : 

» 
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(^rP) t ; t p 2 j^ Y o- Y o2 K '>[ 


Y 02 k r'(02) + Y 0 k r(0) 


+ <*N2 + Y NO + Y N > k '(N2) ] ’ < Y C2 Y N ’ Y NO Y 0 > k Y 


+ Ay 2 _ v Y K^) 

^NO Y 02 Y N2 c ' r 


(2.57) 


(”^)j p2 { (pY N - Yn 2 K c I) [ 


V l(^ + Y 

Y N2 r(N2) Y N r(N) 


+ ^02 + Y NO + Y 0 ) k r (02)] + ^NO Y N " Y N2 Y 0 K c * k r 


+ (Y NO ” Y 02 Y N2 


K. V ' >C } 

»’{ w 0 y <r '«o r -'" , ' , o2 Y ,» 


(2.58) 


+ ^02 Y N- Y NO Y 0 K c V)k I V - ^NO Y N- Y N2 Y 0 K c’ > k r' 


+Y NO + Y 0 +Y N^ k r 


V V 


III 


" 2 ^NO " Y 02 Y N2 


*»') 


(2.59) 


Proper rearrangement 1 of terms in the corresponding expressions for the mon - 


atomic species O and N shows that 



These equations con be obtained, in a simpler way, directly from the moss balance 
equations, Eq. (2. 19) and (2.20). 

It is Important to keep in mind that the expressions given by Eq. (2.57) 
through (2.61) are equations representing rates of production due to chemical 

temperature . As pointed out by Vincenti and Kruger (Ref. 24), they ore customarily 

extended to open systems of varyir* volume and temperature by identifying density 

and temperature os the instantaneous values for the system. This implies that the 

internal rate processes in a moving fluid are the same as those which occur at the 

some state in a static system. The (d X. /dt)^ is related to the mass rate of 

formation of species i, as defined in the species continuity equation, by the 
relation 


/ d X. v 

w i = M i (~af~) 

' 'tot 


(2.62) 


2^4.4 Reaction Rate Constants 

In order to evaluate the reaction rates given in the previous section, 
expressions are needed for the rate constants. Although theories exist for the 
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calculation of dissociation and recombination rate constants (Ref* 24 , 26)# they 

are either too complex to be useful in practice, or too simplified to produce 

acceptable agreement with experimental data. However, it was recognized on 

experimental grounds that, for many reactions, a plot of the logarithm of the 

dissociation rate constant versus reciprocal temperature produces an essentially 

straight line of negative slope, which leads to the Arrhenius equation 

_ J\_ 

k d = A e k T (2.63) 

where e. is the activation energy and A some constant independent of tempera- 

ture. A refinement of this expression is usually obtained by assuming 

0 * 

k , - A T e T (2.64) 

d 

where 0 * is the characteristic temperature of dissociation, and the constants A 
and n are determined by curve fits to the experimental data. Expressions of this 
type are numerous in the literature and can usually be traced to the papers by 
Wray (Ref. 25) and Lin and Teare (Ref. 26). Expressions used in the present calcu- 
lations were taken from Spurk, Gerber and Sedney (Ref. 20), and are listed in 
Table III. 


2.5 Equilibrium Constants 

The equilibrium constants, defined in Eq. (2.47), are also needed for the 
evaluation of the reaction rates* Following again Vincenti and Kruger (Ref. 24), 
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It is found that, in terms of the partition functions, the equilibrium composition of 
a gas is given by 


v. 


TT n. ' 


TT v * 

1 1 n. i 
i • 


n 


v ' 


[ . L Q 


v. 

TJq ' 


exp 


AE 


* 

R T 


(2.65) 


where AE is the difference in the bond energy of the reactants in Eq. (2.40), n. 
are the numbers of particles involved, and v. and X' are the stoi chic metric 

i i 

coefficients. The n. are easily expressed in terms of the concentrations by 
noticing that, according to Eq. (2.41), 


C. 

X = L 

i M. 

i 


P. 


P = 


M. 


n . 

i 

VNT 


( 2 . 66 ) 


Therefore, expressing the left hand side of Eq. (2.65) in terms of the concentrations 


X. , it can be shown that 

i 


V.' - yv v' 

X.’ 1 1 TTn. ' 

1 1 = (N*V) 1 


v. 

nx . 1 

i » 


v. 

nn, • 


(2-67) 


Combining this with Eq. (2.47) and (2.65) it is seen that 

TT q v ! 

K c (T) = (N*V) — ! ^ 

TJ Q 1 


exp 


A E 
R* T 


( 2 . 68 ) 


which, when applied to the reactions given in Eq* (2*51) through (2.56)/ turns out 
to be a function of temperature only. Eq. (2*68) shows also why the contributions 
due to electronic excitation must not be neglected in the equilibrium constants. 
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Here# the ratio of the partition function is important, while in the internal energy 
their contribution is only additive. 

It is clear that substitution of the proper expressions for the partition 
functions in Eq. (2*68) leads to complicated functions of the temperature. In 
particular, since use is often made of the equilibrium constants in comparing 
dissociation and recombination rate constants, it is convenient to have them in a 
simple analytic form. For this purpose Wray (Ref. 25) has fitted the equilibrium 
constants calculated from Eq. (2.68) for the reactions given by Eq. (2.51) through 
(2*56) with Arrhenius type expressions of the form 

K (T) = A T " exp (- 9*/T) (2.69) 

c 

These expressions, given in Table IV/ are employed in the present calculations. 

With respect to these curve fits a remark is in order at this time. As men- 
tioned earlier, the principle of detailed balancing requires that, for a gas mixture 
to be in overall equilibrium, every molecular process and its inverse must be 
individually in balance. Thus it can be shown that, if the total reaction rates and 
their individual contributions are set to zero, the concentration equilibrium con- 
stants, for the reactions chosen herein, result as follows: 


K 

K 



c 

II _ 


c 


< X 0 / X 02 } equil. 
(X N^ X N2 ) equil. 


(2.70) 


(2.71) 
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K " U < X 0 X N /X N0>e q uil. 


and It can also be shown that 


IV _ 


K 


III 


(2.72) 


(2.73) 


V _ 


K 


II 


K 


III 


c 

v, K ' X” 

„ VI _ c c 

r\ 1 " 

c in 2 

K ) 


K 


IV 


VI 


(2.74) 


(2.75) 


Since Wray curve fitted all the K except for K , which was determined from 

c c 

K ^ and according to Eq. (2*75), it is not surprising that Wray's expressions 
c c 

for and differ from and as given in Table IV. For the 

c c c c 

I III 

sake of consistency, the latter expressions were obtained from K through 

according to Eq. (2.73) through (2.75). Numerical calculations show that the 

disagreement between the K given by Wray and those obtained from equations 

c 

(2.73) through (2.75) is negligible. 


2.6 Equilibrium Composition 

It was mentioned at the beginning of this chapter that is is of interest to 
calculate equilibrium flow even though it is, in the exact sense, a purely 
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hypothetical case. In order to do this, the equilibrium composition of the mixture 
has to be determined as a function of other flow variables. 

For a five component mixture this means that five linearly independent 
equations are needed to solve for all the unknowns. Equations (2.18) and (2.19), 
together with Eq. (2.70) through (2.72), provide just such a set. Hence, from 
Eq. (2.18) and (2.19) : 


L C. = 1 


(2.76) 


M 


CU +C 


O 02 \ M 

while from Eq. (2.70) through (2*72): 


(m wi ^ ) 


NO 


2 P C, 


NO 


- m 


02 


M n K 
O c 


I 


(2.77) 


(2.78) 


2 P C 


N 


N2 K 11 

N c 


(2.79) 


NO 


pCq c 


N ( 


M 


NO 


M o M N 


) 


(2.80) 


Substituting the expressions for the mass fractions Cq 2 > an ^ into the 

mass balance equations yields 
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q + q + r ^ ■+ . Z £ r 2 

° N M n K> ° M k ,K n N 

O c N c 


and 


('- NO ,., \ c c 

\ M_ M. , / '-O N 


K 1,1 V M C M N 
c 


= 1 


Q + r ^ 

° M 0 K' ° M M K m "°" N 

O c N c 


c 


m 


from which C.. is obtained In terms of C~ as 
N O 


M r 111 

„ _ m n k c r 

N ?C n L 


K ' - ^ 20 ,-2 
•o L m ° °'m„k' ° 


] 


Using this equation to eliminate from Eq. (2.81) eventually results i 


fourth order equation for C^: 


b 4 C 0 + b 3 C 0 + b 2 C 0 + b 1 C 0 +b 0 0 


b n = 2 (in M_ ) 2 K_ 1 K UI 


c c 


b , ~ m M 
I o 




(2-81) 


(2.82) 


(2.83) 


(2.84) 


where 
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b 2 =M 0 


1 -m 


b 3 =PM 


Q K 1 K V [ — 

O c c l k 


O c 


“ 2 K Iv - ' 
Vi c 


M 


N 


] 


b ., ( !k V [- 

4 c L K 


VI 


- 1 


] 



(2.85) 


Since the b's are functions of temperature and density only, Eq. (2.84) determines 
the equilibrium composition for given temperature and density. It is seen from 
Eq. (2.77) that the only possible solution can be selected by requiring the atomic 
oxygen mass fraction Cq found from Eq. (2.84) to satisfy 

0 SC 0 S [-o- C 02-(^) C Ho] < 2 - 86 > 


2.7 Review of Assumptions 

It is instructive to examine the validity of the assumption of vibrational 
equilibrium which was employed throughout this chapter. A measure for the 
validity of this assumption is cpven by the ratio of the chemical over the vibrational 
relaxation time of the respective species. 
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Using expressions given in Ref. 27, vibrational relaxation times for 
and Nj are given by 


v 

02 


t_ p = 1 .6212 * 10 


-4 


exp 


( 101.44 \ 

\ T 173 “ ) 


(2.8 7) 


T N2 P = 1 * 1146 * 10 " 6 j] ^ 2 ex P 


(>) 


( 2 . 88 ) 


where is given in sec for the pressure in N/m and the temperature in K 

« 

As previously indicated/ the quantity ( p/W. ) has the dimension of time. 
A chemical relaxation time can thus be defined as 



where the mass rate of formation/ 


P 

* 

W. 

i 

W. / can be evaluated for and Nj 


(2.89) 


from 


Eq. (2.57), (2.58) and (2.62). 


Ratios (t*V t v ). for i = O ^ are shown in Fig. 2 as function of 

temperature, together with results for O ^ taken directly from Ref. (28). It 
should be noted that the values of T d / T V are independent of density. 

It is seen from Fig. 2 that, at least for molecular oxygen and nitrogen, 
vibrational relaxation times are always shorter than chemical relaxation times for 
temperatures below 12,000 °K. Hence it is concluded that, for the temperature of 
interest in the present investigation, vibrational relaxation will in general proceed 
much faster than chemical relaxation. Therefore, the assumption of vibrational 
equilibrium is a justifiable simplification. 
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However, it is also recognized, that for temperatures above 8,000 °K 
vibrational and chemical relaxation times assume values of the same order of 
magnitude. Strictly speaking, vibrational relaxation and vibration-dissociation 
coupling should therefore be considered in this temperature range. The difficulty 
here is that, although vibrational relaxation rates are known, an accurate model 
for the coupling does not seem to be available. Conclusions drawn by investi- 
gators in this particular field. Ref. 19 for example, leave the question open as to 
the relative validity of the various models. Rather than to introduce equations of 
which the effect is not yet well understood, the vibrational excitation is assumed 
to be uncoupled to the dissociation, and in equilibrium at all times. 



CHAPTER III 


FORMULATION OF THE PROBLEM 

3.1 Basic Equations 

Neglecting body forces, viscosity, diffusion, heat conduction, and radi- 
ation, the basic equations for sreudy, adiabatic continuum flow are the following: 
Conservation of mass: 

7 • (p q ) = 0 (3- 1) 

Conservation of momentum: 

2 - _ i 

7 (_a~) +(7 X q)xq+j-Vp-0 (3.2) 

Conservation of energy: 

2 

q • V (h + -?j— ) ~ 0 

or, assuming all streamlines to originate from the same reservoir, 

2 

h + = H = constant (3-3) 

These equations are the common equations of motion which do not depend on any 
particular gas model • In order to solve these equations for the unknowns, they 
must be supplemented by a thermal and a caloric equation of state which do depend 
on the particular gas model and were derived in Chapter II as 

p = pRZT (3-4) 
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and 

h = E C. h. fr) (3.5) 

i 1 ' 

If the gGs under consideration is chemically reacting, additional unknowns arise 
in form of the mass fractions, C. , hence, additional equations for the conser- 
vation of species must be added to the system. For the air model used herein 
these are given by Eq. (2.39), which, for steady f low, reduces to 

W i 

q • V C. = — (3.6) 

^ i P 

where i = O, N, and NO, and the mass balance equations given by equations 
(2. 18) and (2.19). 

The choice of the coordinate system is largely a matter of experience. For 
the problem at hand, an orthogonal curvilinear coordinate system with coordinates 
tangential and normal to the body surface (see Fig. 3) seems to be the most suitable. 
From Ref. 2° one obtains for the gradient 


VF = — If- + ° 2 5F + li JLL 

h, iq Sx 2 h 3 3 x 3 


(3.7) 


where F is any scalar point function. For a vector point function F the diver- 


gence is given by 


h 1 h 2 h 3 


2 3 (3.8) 


and the curl by 
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Vx F = 


1 

h 1 h 2 h 3 


The metric coefficients h ^ 


s ; > 

h 2 e 2 

h 3®3 

d 

d 

d 

3x, 

dx 

d x _ 

2 

3 

h ,F, 

h 2 F 2 

h 3 F 3 


h^ follow from the differentia! arc length 


(3.9) 


(ds) 2 - (h^ dx ^ 2 + (h 2 dx 2 ) 2 + (h 3 dx 3 ) 2 (3.10) 

Using as coordinates x, y, and (1), where <t> is the circumferential coordinate# 
one obtains from Fig. 4 

h j dXj = ( 1 + Ky ) dx 

h 2 dx 2 88 d y 

h 3 dx 3 - r^ d <t> = (r^ + y cos 0)^ d <t> (3*11) 

where K = l/R is the body surface curvature, and 

{ 0 for plane flow 

1 for axi symmetric flow 

Applying equations (3.7) through (3. 1 1) to equations (3.1), (3.2)and (3.6)and 
restricting to plane or axi symmetric flow, that is, zero angle of attack, the 
resulting equations are as follows. 

Conservation of mass: 


■^-(pur i ) + -^[(1+Ky) Pvr j ] = 0 


(3.12) 
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x - Momentum: 

u -AiL + ( 1 + Ky) v -1^ + uvK + — 4*-- 0 (3.13) 

dx ' oy p ox 

y - Momentum: 

u |^ + (' + K y )v-^- - u 2 k+ -ii&- -§■£-= 0 (3.14) 

Conservation of species: 

dC. SC. W. 

“•sT 1 +(,+Ky)v T7 ‘ (,+Ky) -F" = 0 (3 - ,5) 

The above equations/ together with the thermal equation of state, the energy 
equation, and two mass balance equations, form a set of ten equations for the 
unknowns u, v, p, p, T and five C.. 

Before discussing the details of the solution of the above system, it will be 
convenient to have equations (3. 13) through (3. 15) also available in the divergence 
form* After some algebraic manipulations the following equations can be obtained, 
x - Momentum: 

^-[(p + Pu 2 )r j ] + -g|-[(l+Ky)puvr J ] 

+ p u v K r 3 - j p ( 1 - Ky) sin 9 - 0 (3.16) 

y - Momentum: 

-JL- (puvrV*^-[ (1 + Ky) (p + p v 2 )r 3 J 

- (p + p u 2) K r 3 - j p ( 1 + Ky) cos 9 52 0 (3*17) 
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Conservation of Species: 

-A- (puC.iVf^[ (l+Ky)pvC 1 r J ]-(l + Ky)r 1 W ( = 0 (3.18) 

The boundary conditions which are needed to solve the equations established here 
are given in the next section. 


3.2 Boundary Conditions 

3.2.1 Body Surface and Shock Conditions 

The body surface is assumed to be chemically inert . The condition for flow 


tangency on the body surface is given by 

v b = 0 (3.19) 

The conditions behind the shock wave are obtained from the conservation of mass, 
momentum and energy across the shock wave. 

Conservation of mass: 



V 

n s 


(3.20) 


Conservation of momentum: 


V = V 
t 00 ts 


^ n 00 ‘ Poo 


P V + P 
s ns s 


Conservation of energy: 




00 


= h + v 2 

s 2 ns 


(3.21) 

(3.22) 


00 


(3.23) 
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where the pressure and the enthalpy are given by the equation of state, Eq. (2. 16), 
and the enthalpy equation, Eq. (2.30), respectively. It should be noted that the 
form of the enthalpy equation implies that the molecular vibrations are in equili- 
brium everywhere. 

For the calculations to follow, the velocities behind the shock must be 
given in the present coordinate system. The following geometric relations are 
obtained from Fig. 5: 

u = V sin 3 + V. cos 3 (3*24) 

s ns t s 

v - - V cos 3 + V. sin 3 (3.25) 

s ns ts 

where the velocity components as referred to the shock are given by 

p oo 

V = V_ sin cr (3.26) 

ns p 00 
s 

= V t „ = cos a (3.27) 

Fig. 4 cl so indicates that the relation between the shock coordinates and the 
shock wave angle cr is given by 

d y 

■■■ " 55 ( 1 + Ky ) tan 3 (3.28) 

dx 's' 

where, in general, 3 is a function of x. 
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3.2*2 Frozen Shock Conditions 


For frozen and nonequilibrium flow it will be assumed that the composition 
of the air does not change across the shock# that is# Z B “ Z^_. Since# in the 
present calculations# the shock wave is assumed to be a discontinuity of zero 
thickness# and any possible chemical reaction would need some finite thickness 
in order to take place, the assumption is a logical one* 

Combining the equation of state with equations (3*20) through (3.23)# it 
can be shown that 


T 

s 


n “ 
R Z 



2(h s -KJ- 


n‘ 


-2<h s -hJ 


R Z 


(3.29) 


For any given free stream conditions (i *e* T^# V ro # C.^) and the shock wave 
angle c# the temperature behind the shock# , can now be found from equations 
(3*29) and (2.30) by solving these equations numerically. Using the continuity 
equation# Eq. (3.20)# and the energy equation, Eq. (3.23), the density ratio 
across the shock is obtained as 



Knowing the temperature and the density 
ratio across the shock becomes 



<h s -hj 


(3.30) 


behind the shock, the static pressure 



(3.31) 
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which completes the calculations of the static conditions on the downstream side 
of the shock wave . 

Unfortunately, the solution described above involves an iterative calcula- 
tion of the temperature T s . Since iterations are normally time consuming, it is 
impractical to calculate the shock conditions by this process at each step in the 
course of a numerical calculation. Instead, the expressions for u^ , v^, T^, , 

and p , presented in this and in the previous section, may be expressed in diffe- 
rential form in terms of the shock wave angle gradient. It is then possible to 
integrate for some selected key variables, say density P g and temperature T s , 
while the values of the remaining variables are determined from simple algebraic 


equations < 


Differentiating the expressions for u^, v^, T^, P^ and p^ with respect to 


x results in the following differential shock relotions: 


d T J rr 

- = T n, -£E- 

d x 00 1 d x 


( 3 . 32 ) 


dp 


$ „ 


d x 


n da 

°2 dx 


( 3 . 33 ) 


d x 



( 3 . 34 ) 


d u 
s 

d x 


= v„n 3 


d a „ 

d x s 


( 3 . 35 ) 
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d v 

L = V 0 

d x V “ 4 d x 


d a 


- u K 
s 


(3.36) 


where 0 through 0 are dimensionless functions explicitly given in Appendix 
1 4 


B. 


From Fig. 3 it may be deduced that 

r 55 r, + y cos 9 
s b 's 

which, upon differentiation, yields 


d r 


s _ 


" ( 1 “ Ky ) sin 0 + ( 1 + Ky ) cos 9 tan 0 
dx s' 's' 


where use was made of the relation 


d r, 


d x 


b _ . a 
- sin 9 


(3.37) 


(3.38) 


(3.39) 


which is obtainable from geometry 


3.2.3 Equilibrium Shock Conditions 

In this case it will be assumed that the gas instantaneously reaches the state 
of thermodynamic, that is, thermal and chemical, equilibrium behind the shock, 
irrespective of the conditions in the free stream which may be arbitrary. While 
for frozen shock conditions the composition of the gas on the downstream side of 
the shock wave is prescribed by the free stream composition, it is unknown in the 
present problem. The necessary additional relations are furnished by the equations 
for chemical equilibrium, derived in section 2.6. 
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Using the equation of state to eliminate the pressures from the equations of 
conservation of mass and momentum across the shock/ it can be shown that 


where 


00 


2 




RZ T 
s s 


(3.40) 


n 


00 


b = 1 + 


RZ^T* 


(3.41) 


V 


n 


The reciprocal of Eq. (3.30)/ which was derived from the conservation of energy. 


yields 


f- 


2 (h -h ) 


(3.42) 


00 


For given shock wave angle c, equations (3.40) and (3.42) are two independent 

relations for the density ratio across the shock as function of density and 

temperature T . Hence, assuming some reasonable values for p and T , a 
s s s 

solution is obtained (by iteration on p and T ) when the difference between the 

s 5 

right hand sides of equations (3.40) and (3.42) assumes a specified minimum value. 

Having determined the density p and the temperature T , in the course 

s s 

of which the species mass fractions C. have also been calculated, all other 

I 5 

flow variables are readily evaluated from the available equations. 

Although equilibrium shock data are available in the literature (Ref. 30, 31), 
these data are usually presented for selected atmospheric free stream conditions at 



42 


relatively low altitudes. The procedure described above, however, is applicable 
for arbitrary combinations of free stream pressure, temperature and composition, 
as they may occur at high altitudes or in a hypersonic wind tunnel . 



CHAPTER IV 


METHOD OF SOLUTION 


4.1 The Method of Integral Relations 
4.1.1 The Method in General 

In order to apply the method of integral relations as described by Dorodnitsyn 
in Ref. 2, all the partial differential equations are cast into the divergence form, 
namel y 


3 F. a G . 

L_ + i 

dx by 



(i * 1, 2, m) 


( 4 . 1 ) 


Here x and y are the independent variables, while F. , G. and H. are the 

i i i 

known functions of the dependent variables, and m denotes the number of 
equations. 

Consider now the solution of the above system of m equations in the region 

which is bounded by the body surface and by the shock wave y (x). Dividing 

s 

this region into N curvilinear strips (see Fig. 6) bounded by lines 

v = y k = -jq- y s to (4-2) 

where k = 0, 1,2, . . ., N, the system (4.1) can be integrated with respect to 
y across each strip. The result is a system of integral equations of the form 
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y k + l aF y k + l 

f -W- dy+G i,k + l- G i,k + / H i dy = 0 (4.3) 

y k y k 


where k - 0/ 1, 2, . N * 1. Using Leibniz's rule for differentiation under 
the integral sign (Ref. 32), the first integral may be rewritten and the above 
equations become 



+ G. 


i,k + 1 



y k + l 

+ f H. dy = 0 
y k 


(4.4) 


In order to evaluate the integrals, the integrands must be known functions of y. 

Unfortunately the integrands F. and H. contain precisely those dependent 

variables for which values are to be calculated, hence some approximations 'iave 

to be made. Generally, as an approximation, any interpolation formula can be 

used which expresses the value of F. or H. at an arbitrary O^y 221 y through its 

II s 

values on the lines y = y^ . For example, one may use polynomials of the form 



(4.5) 
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where the a. (x) and b. (x) are to be evaluated in terms of the boundary values 
in in 

at the strip interfaces, say 

a ty=y k : f. = F. k 

H. - H |ffc (4.6) 

Performing the necessary operations on these polynomial approximations yields a 
system of m. N ordinary differential equations of the following form: 


N 



n -0 


where k = 0, 1 , N - 1; n = 0, 1, . . . , N and i - 1 , 2, . • . , m. 

The solution of the system (4.7), obtained when the region is divided into 
N strips, is called the Nth approximation. It is obvious from the above that 
with higher approximations the complexity of the problem increases considerably . 
Belorserkovskii (Ref. 37) demonstrated the convergence of the method by comparing 
one**, two-, and three-strip solutions for supersonic flow past a circular cylinder. 
He showed that, at least for that case, there is practically no difference between 
the three solutions. Moreover, since previous investigations with simple rate 
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processes (Ref . 12, 13) have shown that the first (linear) approximation also agrees 
well with results from other methods, the linear approximation is investigated 

herein. 


4.1.2 The First Approximation 

For the first approximation, according to the previous section, N 85 1, k“0, 
n * 0, 1 . Hence from Eq. (4.5) the polynomial approximations have the general 

form 

P. - c. 0 (x) + c. j (x) y (4.8) 

According to (4.6), and the coordinate system, (Fig. 3), it is found that at 

y = 0 ; P i,b ” C i 0 ^ 

(4.9) 

y s V P i, s = c io w+c n wy, 


From (4.9) the coefficient functions in general are 

c i0 (x) = P i/b 

•llW = T" (p !,s- p i,b) 

$ 


(4.10) 


Applying this procedure to F. and H., the coefficient functions are given by 
Eq. (4.10). Equation (4.7) may now be simplified and, for the first approximation. 


results ir 
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d 

d x 







+ -2— (G. -G. J + (N. +H. .) = 0 (4.11) 

y i/S i/b i/S i/b 
s 

The obove equation might be called an operator equation, permitting the 
transformation of the given partial differential equations of the form (4.1) directly 
into ordinary differential equations for the dependent variables along the body 
surface and the shock wave boundary. 


4.2 Nonequilibrium Flow 

4.2.1 General Remarks 

As early as 1929, Busemann (Ref. 33) showed that supersonic flow of an 
inviscid gas about circular cones is characterized by the fact that the pressure and 
the velocity vector are constant on coaxial conical surfaces having the same vertex 
as the conical body. A consequence of this is that the shock wave, forming one 
boundary of this flow field, must itself be a conical surface. It will be shown in 
the discussion of the results, that, in the case of nonequilibrium, chemically 
relaxing flow in this case, the time dependency of the chemical reactions destroy 
this self-similarity mentioned by Busemann, even though the flow is considered 
to be inviscid. Therefore, no definite statement can be made, a priori, about the 
shock shape or the behavior of other flow parameters. Values for all variables 
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must be calculated by integrating the conservation equations together with the 
equation of state and the boundary conditions in a step-by-step fashion along the 
length of the cone. 

Two different sets of equations will be examined. The first set, later 
referred to as the standard approximation, contains all the conservation equations, 
equations (3 . 12) through (3.15), in approximate form. The second set, later 
referred to as the semi -exact procedure, uses only the continuity and the y- 
momentum equations in their approximate form, while the x-momentum and species 
continuity equations are used in the exact form. 

For a supersonic flow field, which will be considered here, the equations 
are of quasi hyperbolic character and thus form an initial value problem. Initial 
values are given by the frozen flow solution to be discussed in section 4.3, where 
it will also be shown that initial gradients can be derived as functions of the 
initial values. 

4.2.2 The Standard Approximation 

The application of Eq. (4.11) to the conservation equations, Eq. (3.12) and 
equations (3. 16) through (3.18), results in the following ordinary differential 
equations where use was made of the fact that v^ ~ 0, and y _ 0, at the body 


surface . 
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Continuity: 


( Pl u l r J) + "T"“ (P u fJ ) 
dx bbb dx sss 


-f (P, u , 'l - p b °b r b > TT - j- < ' 1 +K ^> P s v s '! 
$ $ 


x-Momentum: 


ir[ ( Pb + l u b 2)r b] + -h [ (p s +p s u s 2)r s J ] 


r[ (p s +p s u * 2)r s J * <p b H u b 2)r b ] 


dy 5 

dx 


— ( 1 + Ky ) P u v r 1 - p u v K r J 
y ' 's s s s s sss s 
s 


+ j sin 9 [p $ (i -Ky s )+ p^ ] 


y-Momentum: 


(P. u v. r • ) 


J* = 


d vr v ? i 

X s s s s 


i * dy 

— (P u v r J ) 
y ' s s s s dx 
's 


[(! + ■*,> (p.+P.v, 2 ) r> - p b rJ ] 


+ K [ (p s + P s u s 2)r s + (p b + i u b 2)r b J ] 
+ j cos e [ P s ( 1 + Ky s ) + ] 


(4.12) 


(4.13) 


(4.14) 
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Species continuity: 


dT- (p b u b C ib r b J)+ dT (p s °s c i* r * J) 

i • . dy 

= (p y Q m r - * - a u. C . i r, ^ ) — r~ 

y ' s s i s s T) b ib b / dx 


C l. r . ,+ < , + ^.> r . J ^. + r b *ib < 4 - ,5 > 


Expanding the differentials and introducing the differential shock relations 

mm 

given in section 3.22, the above equations assume the following form, where 6 


is defined as 


Continuity: 


6 = 


du, 


y, to 


(4.16) 


ul + Pi -t-”- + ( 1 + j $ cos 0) (p u + P V„ 

b dx b dx J /x * s 2 s “ 3' dx 


_ 1 


^sV^bVVsV 


(4.17) 


x -Momentum: 


^ V ^ + 2 P b u b IST +( 1 ^ 9) [ R Z s < P s T - n 

+ (V.T s n 2 ) + P„u 5 2 ^+2P s V„ Uj C2 3 ] 

{ <Ps + P s U S 2)F I- ( Pb +P b°b 2 ) F 2 

-p u v ( 1 + j &cos 0) ( 2 + 5 Ky ) 

$ s s s 

+ [ P b +P s ^ +Ky j] J ^ s?n9 | 


6 % 


U b 


d a 
dx 


= J_ 

y. 


(4.18) 


y-Momentum: 


(1+j 6cos0 ) (p* u $ v s V P s V « v s ^ +p s V « u s ) ~Tx 
“ -7-| p s v s (u s F r v S F 3 ) + (p b' p , )(2+j fco$ 9+Ky s ) 


+ [ p b u b^ + 2 P s ( 1 + j * cos 0) J Ky s | 

(4.19) 

Species continuity: 

dG. 

P b u b dx " (C ib'Ss )(,+J 5cosS)(p “ U s n 2 


+ p V Q, ) -42- 

M s 00 3 ; d X 

s 

* W. b + ( 1 + Ky $ ) (1+j 6cos 0) W. $ 


L. (c. - C ) (p o F.-p v F*) 

y ' ib is s s i s s 3 

s 

(4.20) 

where F., F 0 , and F^, are dimensionless functions defined as 
12 3 


F j - ( 1 + Ky^) tan 0 - ( 1 - Ky^) j &sin 0 

(4-21) 

Fj ~ (1 + Ky^tan 3 + j 6 sin 8 

(4.22) 

F_ = (2 + 3 Ky )(1 +j 6 cos 0) 
3 s 

(4.23) 


The equations given above must be supplemented by the energy equation and the 
equation of state in differential form. 


Energy: 
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du b 


°b d* 


+ c 


-' T b 

pb d x 


E-, 


d C. 


<b 


b d x 


= 0 


(4.24) 


State: 


dp b 

sr 


- R 


b d x 


- R 


, ^ 

b dx 


- P b R * T b 


Y— 

Lu M. 


d C. 


ib 


d x 


= 0 


(4.25) 


Since the mass fractions C. ere linearly related to each other by the two mass 
balance equations/ only three of the five C. are independent. Therefore/ 
equations (4.17) through (4.25) constitute a system of eight simultaneous diffe- 
rential equations for the gradients of p^ , u^, p^/ C7 , and three of the five 



It may be seen from the conservation equations in the divergence form/ and 
from the discussion of the integral method/ that the equations presented in this 
section contain linear approximations for the following terms: 

Continuity: pu r ^ 

2 j 

x-Momentum: ( p + p u ) r , 

puvKr^- j p ( 1 - Ky) sin 9 
« 

y-Momentum: p u v r ^ , 


(p+ pu 2 ) K r ^ + j p ( 1 + Ky) cos 9 

Species continuity: 

(pu C. r J ), (1 + Ky) W. r J 
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It is of interest to examine these approximations in some detail . First, it is 

noted that it has no effect on the x-Momentum equation, Eq» (4.13), whether one 

approximates the terms p u v K r * and j p ( 1 - Ky) sin 6 separately, or combined. 

The situation is analogous for the y-Momentum equation. Second, the body sur- 

face curvature K (x), and the surface inclination angle 0(x) are constants as 

far as the linear approximations in y-direction are concerned. It is then noted 

v 2 j 

that both of the momentum equations contain linear approximations for ( p + p u )H 
and p u v r^ . However, while the term p ( 1 - Ky) is linearized in the x-Momen- 
turn equation, it is p (1 + Ky) in the y-Momentum equation. This is an obvious 
discrepancy which requires particular attention. 

Application of equations (4.8) through (4. 10) to these terms leads to the 
following relations for p (y) : !' 


p, M - 



p 5 ( 1 - Ky s ) 

k y s 



p 2 (y) 35 



t 'p,(l + Ky,)-p b 
y s 



Since two different values for the pressure at any one point in the flow field are 
unacceptable, p^ (y) and p ^ (y) must be required to be the same. It can be seen 
at once that this condition is satisfied for K - 0. On the other hand, for K 4 0, 
it can be shown that this not only requires p^ and p^ to be the same, but also 


that 
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P (y) ~ ~ P $ ~ constant 

Fortunately this discrepancy is the only one, since the application of 

equations (4.8) through (4.10) to the remaining terms generally leads to only one 

2 

nonlinear function of y, each, for pu , pu v, p + pu , pu C., and W.. A 
procedure to avoid the discrepancy altogether, and simultaneously reduce the 
number of functions for which approximations have to be introduced, is discussed 
in the next section. 

Apart from this it should be noted that certain terms disappear in the con- 
servation equations when either j or K, or both, are zero. In this case, no 
linearizations are introduced for the respective terms. Thus, for example, when 
the equations are applied to wedge flow (K - 0, j - 0), the linearized terms 
which are actually used reduce to (p u), p, (p u ), (p u v ), (p u C.) , and W. . 

4.2.3 The Semi -exact Procedure 

The general form of the conservation equations for the one-strip approxi- 
mation, Eq. (4.11), already indicated that the equations are formulated in terms 
of the variables along the shock wave and along the body surface, only. This 
situation suggests the use of the flow tangency condition, Eq. (3.19), in order to 
simplify some of the partial differential equations before they are integrated in 
y-direction, and thereby converted to ordinary differential equations. As a 
matter of fact, it is seen from equations (3.12) through (3.15) that the x-Momentum 
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equation and the species continuity equation reduce to ordinary differential 
equations in terms of the variables along the body surface when v is set to zero. 
Therefore# only the overall continuity equation and the y-Momentum equation 
need to be used in approximate form, which results in a much simpler set of 
equations for the semi-exact procedure. This set then consists of equations (4. 17), 
(4.19), and (4.21) through (4.25), while the x-Momentum and the species con- 



(4.26) 

(4-27) 


respectively. Thus, the introduction of the above exact equations eliminates the 

• • • 

need to approximate the terms pu C. r J , ( 1 + Ky) W . r ■* , and p ( 1 - Ky). In 
particular the discrepancy discussed in the previous section is removed by not 
using an approximation for p (1 - Ky). 

A careful comparison of equations (4.26) and (4.27) with their approximate 
counterparts reveals that they may be expressed in a combined form by introducing 
certain indicators, say e and cl. Setting the indicators to zero results in the 
exact equations, while assigning the value unity yields the approximate equations. 


Hence, 
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x-Momentum: 


d Pk 9 d Pk du k 

inr +eu b TT + < , + e > p b U b-1T 

+ e(1+j &cos9)[kZ (p T„ 0 , + P„ T s 0 2 ) 


+U S (P.« S n 2 +2P $ V “ n 3 ) ] 


d o 

d x 


(4.28) 


m e 

y. 


j(p s + P s “ s 2 )F,-(P| J +P bUb 2 )F 2 


- p u v ( 1 +j 6 cos 0 ) (2 + 5 Ky ) 
$ s s J s 

+ [pb + P s ( 1 + Ky s ) ] j 6 s in 8 | 


Species continuity: 
dC.i 


p b °b -dT ’ a(c ib ' Ss) < ’ + J‘ J cos 9 Hp- “ s ^ + p 5 n 3 ) -If 


= W + a[( 1 + K y s ) ( 1 + j 6 cos 6 ) W. $ 


(4.29) 




(C ib ~ C ! S ) (p s°s F I " P s v s F 3>] 


Equations (4.17), (4.19), (4.24), (4.25), (4.28), and (4.29) now represent 
a system of six simultaneous equations for the gradients of u^, f^, p^, T^, C.^ 
and cj . Noticing that the y-Momentum equation and the species continuity 
equations yield the derivatives of the shock wave angle o and the mass fractions 
Cjb directly, the above set is easily solved for the remaining gradients. Using 
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the energy equation/ Eq. (4.24), and the equation of state, Eq. (4.25), to 
eliminate the derivatives of u^ and p^, the following equations are obtained; 


^b. _ _ 

i 

< 

CM 

CM 

O 

CM 

< 

CM 

O 

(4.30) 

d x 

°11 a 22 

CM 

a 

CM 

0 

1 

dT b . 

°n A 2 

< 

CM 

0 

1 

(4.31) 

d x 

a ll a 22 

' a 12 a 21 

da _ 

A 3 


(4.32) 

d x 

°33 


dC ib = 
d x 

A 4 
a 44 

(i = 0, N, NO) 

(4.33) 


The coefficients a., and the terms A. are functions of the variables along the 

ij i 

body and the shock wave. They are explicitly given in Appendix C in their 
most general form, valid for any combination of the appropriate values for j, K, 
CL , and e . For the present investigation which is primarily concerned with the 
application of the semi-exact procedure (a- 0, e= 0) to the flow about circular 
cones (j “ 1, K - 0), this implies the use of nonlinear approximations for p u and 
pu v , 


Pu 


p b u b r b 


*( p s u 

r b 


r 

s s 



-^» b r b )(y/y 5 ) 

cos 9 


(4.34) 


p u v r 
s s s s 

r^ + y cos 9 


( 



) 


0u v ° 


(4.35) 
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and a linear approximation for p , 

p (y) 58 p b * — — y ( 4 - 36 ) 

Since it can be shown that the right hand sides of equations (4.34) and 
(4.35) are the sums of infinite alternating series in powers of y / y^, where 
0 £ y /y , the error of the method is determined by the linear approximation, 
Eq. (4.36), for the pressure p (y) across the shock layer. 

Having determined the density , temperature T^, shock wave angle a, 
and the mass fractions C* b (i -O, N, NO) by numerical integration of equations 
(4.30) through (4.33), the remaining C ;b (i * O ^ N 2 ) are calculated from the 
mass balance equations. Knowing the temperature T b and all of the mass 
fractions C fb , the static enthalpy can be evaluated, whereupon the velocity u b 
can be computed from the energy equation. Similarly, with density tempera - 
ture T b , and the mass fractions Cj b known, the pressure p b is given by the 
equation of state. 

As far as the variables along the shock wave are concerned, it suffices to 

integrate equations (3.28), (3.32) anJ (3.33) in addition to (4.32). Once the 

shock wave angle a, density p , and temperature T are known, all other 

s s 

variables behind the shock can be computed algebraic!'/ from the expressions given 
in Chapter III. 

The advantage of the semi-exoct procedure over the standard approximation 
appears not only in its relative simplicity, but also in the fact that no 
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approximations are involved which contain the rate equations or the species ma«? 
fractions. This is important because, for chemically relaxing flow, the C. are 
key variables of particular interest. 

4.3 Initial Solution 

4.3.1 Frozen Flow 

In order to start the numerical integration of the equations presented in 
section 4.2, the values of all variables must be known at x ~ 0. The calculation 
of these values will be based on the assumptions that, first, the flow is chemically 
frozen and, second, that, in case of a body with non-zero curvature ( K + 0) 
there is a small region including the tip where the body may be considered to be 
a circular cone or a wedge (K ~ 0). Since the use of frozen shock conditions 
implies that the flow is frozen at the tip, this assumption is not a new one. 

With the molecular vibrations in equilibrium, and the gas composition 
frozen, there exists no mechanism to destroy the self-similar character of the flow 
field, and therefore the classical result of a straight shock wave and constant 
properties along the body streamline may be applied. Mathematically expressed, 
this means that in equations (4.17) through (4.19), the gradients of all flow 
variables along the body and that of the shock wave angle a can be set to zero. 
This implies that also the right hand sides of these equations must vanish. Hence, 
with 6, defined by Eq. (4.16), taking the form 
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y s _ tan 3 

x sin 6 sin $ 


(4.37) 


and fy ? 2 ' F3 r ®ducing to 

F 1 = (1 - j)tan P 

F 2 - ( 1 + j) tan P 

F^ * 2 (1 + j cot 0 tan P) 


(4.38) 

(4.39) 

(4.40) 


the following algebraic equations, valid for chemically frozen flow about wedges 
(j -0) or circular cones (j - 1), are obtained from equations (4.17) through (4.19). 
Continuity: 

P s u s “ i) tan P “ % u b 0 + J*) tan ^ 

-2p v (1 + j cot 0 tan PJ^O (4.41) 

s s 

x**Momentum: 

(P s + p s u s 2 ) (1 " j)tanP- (p^ + u b 2 ) (1 + j ) tan P 
- 2 p $ u $ ( 1 + j cot 0 tan P ) + j (p^ + p $ ) tan P = 0 (4-42) 

y-Momentum: 

2 

P u v (1 - j) tan P - 2 p v ( 1 + j cot 0 tan P ) 
s s s s s 

+ (p^ - p s ) (2 + j cot 0 tan P) - 0 (4.43) 

Equations (4.41) through (4.43), together with the energy equation and the 
equation of state, then represent five equations for five unknowns: u^, p^ , T^p^, 


and o . 
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For a given flight configuration, that is given semi-vertex angle 9 and 
free stream conditions, the calculation is started with an assumed value for the 
shock wave angle cr , and calculating the temperature T , the density p , and 

S 5 

the pressure p from the frozen shock conditions, that is, from equations (3.29) 
through (3.31), respectively. With the velocity components u and v , which 
are to be computed from equations (3-24) and (3.25), the velocity at the body 
surface, u^, is given by 

(4.44) 


u b U s 2 cot (3 + j cot 0 

This relation can be derived by first eliminating (Pj^ " P $ ) from the two momentum 

equations, and then replacing the term (p b u^) from the continuity equation. 

Knowing the velocity u^, the continuity equation yields the density p^. 

With u, and P, available, the surface pressure can be calculated from one of 
b b 

the momentum equations. Finally, the equation of state serves to obtain the 

temperature T^. The energy equation, not needed so far, is now used to check 

the value of u. • If the velocity u, so calculated does not agree with the value 
b b 

obtained from Eq. (4.44), a new value for the shock wave angle o must be chosen; 
and the procedure is repeated. 


4.3.2 Initial Derivatives 

The necessity for determining the derivatives of all variables at x - 0 arises 
from the fact that the equations for the initial values were obtained by equating 
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the numerators of the right hand sides of equations (4. 17) through (4. 19) to zero* 
Since, for on attached shock wove, olso y $ (0) = 0, the right bond sides of 
these equations become indeterminate at the tip of the body. The condition of 
frozen flow (C. = constant) leads to the same difficulty in Eq* (4.20). 

A set of linear equations for the initial derivatives can be derived from the 
governing equations by evaluating their 1 1 mi ti rig expressions as x 0. For this 
purpose let the right hand sides of equations (4.17), (4.28), (4.19), and (4.29) 
be denoted by B y B 2 , B y and B^, respectively. Specializing them for the 
case of zero curvature (K = 0, r^ “x sin 9), and introducing 

6 - 1±- (4.45) 

X 


they may be written as follows: 


B, s 


— Tp u (tan 9- j 6) - a u, (tan 3 +j 6) 

V L S S do 

- 2 P v ( 1 + j 6 cot 0)1 

S S J 


B = 

2 y. 


£ P s u s 2 ( tan 9- j 6) - P^ u b 2 ( tan 0 + j 6) 

- 2 P s v $ u s ( 1 + j Scot 0) + (p s - p^) tan 0 J 


(4.45) 


(4.47) 


B 3 = T 


[p s u s v s (tanP-j5)-2P s v s 2 ( 1 + j 6 cot 0 ) 
-(P s - p^) (2 +j Scot 0)] 


(4.48) 
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B 4 = W. b + a | (1 +j Scot 9) W. s 

L(c -C. ) Tp u (tan 3- j &) - 2 p v ( 1 + &j cot 0 )] 1(4.49) 

ib is L s s s s JJ 


Then applying L ‘Hospital 's rule in order to resolve the indeterminacy, the following 
limiting expressions can be obtained: 


lim B , 

x-q 



P 

s 


V^b 

sin 2 3 


+ 


P v 
s s 


• 2 o 

sin P 


+ 0-J)(P ro « s Q 2 +P s V » Q 3 ) 


“ 2 (cot 3 + j cot 9 ) (p w v $ 0 2 + p $ 

- (,+ j ) T^ (p b u b> (4,50> 


m '{ 


(P s +P s u s 2 ) ~(P b + p b u b 2 ) + P s u 


V 

s s 


sin 2 3 


sin^ 3 


+ (1 -])(P.u 2 °, + 2P V M o 0 ) + RZ (p T„fi, + P„ T s n 2 ) 


2 (cot 0 +j cot 9)(P 00 u s v s n 2 + p s V “ v s n 3 + p s V » u s "d 

- eO + J) (? b “ b 2 )- e “ST 


do_ 

X 


(4.51) 
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lim B 0 

x~*0 J 


{ P U N 
s s 

sin 2@ 


, 2 

P, + P v -ft 
+ —2 LJ r 


sin ^ 0 


+ 0-j) (p„ u, v $ n 2 + P s v„v s n 3 +p s v„u. nj 


$ 4 


- 2(cot0 + j cot0) (p^ v s +2 P s V M ^ 4 )v s 

- R Z $ (2 cot B + j cot 9) ( P $ T. n i + Poo T $ Oj) | 


d a 
dx 


+ ( 2 cot 0 + j cot 0 ) 


"b 


d x 


(4.52) 


and, finally, 

lim B, * W., + a j (1 + j tan 0 cot 9) W. 

x-o 4 ,b L ,s 

"<' + H“b ~TT ] < 4 - 53 > 

where it is noted that the above expressions are to be evaluated at x “ 0 only. 

If the limiting process is now applied to equations (4.17), (4.28), (4. 19), 
and (4.29), and if the expressions given above are inserted on their right hand 
s'des, the following equations for the initial derivatives in terms of the initial 
values result. 


Continuity: 


dp, du, 

u — — + Q + Q v 0 

bdx 1) dx woo 5dx 


da 


* 0 


(4.54) 


x-Momentum: 


Pb 


2 dP b 


0 + *>- 37 - + P( 2 +j) u b "d 


[l + e (3+ 2 j )] P ( 


d u< 


b 1 ^ ~sn 


b + «P_V 2 0. da 


00 00 


6 dx 


= 0 (4.55) 
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/-Momentum: 

4r- + P„ v?, Cl 41 . = 0 (4.54) 

d x 00 co 7 d x 

Speci es conti nui ty : . r 

[)+<X(l + j)]p b u b -gA 

« • 

* W*k + a, ( 1 + j tan 3 cot 9) W. s (4-57) 

where fi 5 ,n and fiy are dimensionless functions listed in Appendix B. 

Supplementing equations (4.54) through (4.57) by the energy equation and 
the equation of state/ 


du. 


P b u b dx 


+ D l!b + p b V -, 

P b C pb d x X b 


0 


(4 58) 


% 

d x 


, dP b dT b 

" RZ b T b ~dT _P b RZ b dx 


p b T b V » 


XT 


00 


S b =0 (4.59) 


where X is a characteristic length/ and L and S are dimensionless functions 
defined in Appendix A, this set can be solved for the initial gradients. The 
resulting expressions are: 


(■&) ; [ ,+e ( ,+ j>] 


1 


u 


00 


XF 


< RZ b L b- 


C pb T b 


S.) (4.60) 


where 
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F= °b 2 { [ 1+S(1 +j) ] e pb-< ,+e > RZ bj n 7 

“ RZ b c pb T b { [' +ep +2 i) ]'X' a 5 - « V 0 +e > n 7 j 
- eRZ b u b [ (2 + j) “5 - n 6 ] 

(-1^) P " V “ -~ ;{ [ ,+e (3 + 2j)]-^- n 

V d x / n +^/i +1)^ 2 ( L J v ® 5 


x “0 [1 +e( 1 +j)]u b 


- c 0 + ( 1 + e) 

6 




(4.61) 


(4.62) 


/ dT b \ . . p - v - ( (\\ h 
' d x i=0 P b C pb I ' P » ' X 


e(2 + j).'Ib. Q -eQ, +(1 +e) CL 
Vce> 5 6 / 

1 +e (1 +j) 




.63) 


/ dC ib> 

W., + a ( 1 + j tan 0 cot 9 ) W . 

t _ ib 15 

V dx J 

* =0 P b « b 

1 + a(l +j) 



(4.64) 


It should be noted that, since no assumptions or approximations were used in 
deriving the initial gradients from the governing equations, the compatibility of 
equations (4.60) through (4.64) with equations (4.30) through (4.33) is assured. 
Again, for a ^ 0 and e= 0, the expressions given above reduce to the initial grad- 
ients to be used for the semi -exact procedure. Letting a - 1 and e " 1 results in 
the initial gradients for the standard approximation- 
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4.4 Equilibrium Flow 

As mentioned earlier, an equilibrium flow solution in connection with non- 
equilibrium flow about wedges or circular cones is useful because it provides a 
check, in an approximate way, for the values which should be approached 
asymptotically by the nonequilibrium solution* There are two reasons for the 
approximation* It will be shown in the discussion of the results that, for non- 
equilibrium flow about cones, the shock wave is curved. By assumption, the 
shock wave angle o' at the apex is the same as that for frozen flow, while further 
downstream, a approaches the equilibrium shock wave angle which is always 
found to be smaller than the shock wave angle a for frozen flow. Consequently, 
those streamlines which have passed through the stronger portion of the shock wave 
form a region close to the body surface which must have a higher entropy than the 
stream lines farther away from the surface. This region is sometimes called an 
entropy layer. Furthermore, it is known from theory that, in contrast to equili- 
brium, nonequilibrium dissociation and recombination are nonisentropic processes. 
Thus, the chemical relaxation in the flow provides another mechanism by which 
entropy is increased* It is therefore concluded that the equilibrium state which is 
reached after complete relaxation has a higher entropy than an equilibrium flow, 
originating from the same free stream conditions, would have. Since, by assumption, 
the two flows have the same total enthalpy, but their entropy is different, they 
cannot reach the same final state- 
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Assuming that the air is in thermodynamic equilibrium everywhere in the 
flow field, the same considerations as for frozen flow apply, namely, that the 
flow field is similar. Again, this implies that all flow field parameters are 
constant on planes, or coaxial cones intersecting at the tip of the body. 

Newman (Ref. 36) has used equations (4.41) through (4.43), in connection 
with thermodynamic data from tables, to calculate conical flow parameters for air 
in equilibrium. He found that these algebraic equations yield results which are in 
excellent agreement with the results from the numerical integration of the Taylor- 
Maccoll differential equation by Romig (Ref- 21). Similar calculations by Thoenes 
(Ref. 16, 17) for a simple air model, but using equations analogous to those pre- 
sented in Chapter II for the thermodynamic properties, confirmed this agreement. 
The more sophisticated air model employed in the present investigation requires a 
method of calculation which differs slightly from that used previously, and is there- 
fore outlined below. 

For an assumed shock wave angle cj the equilibrium conditions behind the 

shock wave can be calculated as described in Section 3.2.3, whereupon the body 

surface velocity, u, , is computed from Eq. (4.44). Substituting this expression 

b 

into the continuity equation, Eq. (4. 41), yields 



[ P s u, ('-]) 

- 2 v (cot 3 + j cot 0 )J 

(2 cot 3 + j cot 0) 

0 +j) 

u s (2 cot 3+ j cot 0 ) + v s 

] 


(4.65) 
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Solving the y-Momentum equation, Eq. (4.43), for the body surface pressure 
results in 


Pb =P s + 


2 p v* ( 1 + j cot 0 tan ^)“P $ U S V S (1 " j ) ^ ari ^ 
2 + j cot 9 tan 3 


(4.66) 


Starting from a first guess for the body surface temperature, , together with the 
density at the body surface, , from Eq. (4.65), the equilibrium species mass 
fractions must now be calculated by an iteration, using th*. equations provided in 
Section 2.6. This iteration can be terminated when the mass fractions C.^ and 
the temperature thus obtained, together with density and pressure p^ 
from equations (4.65) and (4.66), satisfy the equation of state. 

The energy equation, not used so far, again provides the closing link in 
the loop. If the value for the velocity at the body surface, u^ , obtained from 
the energy equation does not agree with the one calculated from Eq. (4.44), the 
procedure has to be repeated by selecting a new value for the shock wave angle cr . 


4.5 Numerical Techniques 

For the numerical evaluation all equations were converted to dimensionless 
form by introducing the following dimensionless variables: 


S 
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u 


f 




9 


P 


i 



h ' 




(4-67) 


Since infinite long cones have no typical dimension/ the characteristic length X 
was defined as a dissociation relaxation length by 


X = ! (4.68) 



where Z is the compressibility factor defined by Eq. (2.17). The equations were 

then programmed in FORTRAN V, suitable for the UNIVAC 1108 located at the 

University of Alabama Research Institute. For simplicity of program check-out 

and debugging the problem was programmed in three separate parts for frozen, 

equilibrium, and nonequilibrium flow. 

The frozen flow program consists essentially of a double loop iteration, 

following the process of solution as described in Section 4.3.1 . The process is 

started with an assumed value for the shock wave angle o in the main program, 

which then calls a subroutine for the iterative calculation of the temperature 

behind the shock, T , from Eq. (3.29). If the difference in the body surface 

s 

velocity, u^ , which is computed from two independent equations, is larger 
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than a specified minimum, the computation returns to the main program and resumes 
with an adjusted value of the shock wave angle <7 . The built-in limit of ten 
iterations was in no case exceeded . Running times, including compilation of the 
program, are in the order of 5 to 10 seconds, even with double precision. 

The structure of the equilibrium flow program is analogous to the one for 
frozen flow, however, with additional loops for the calculation of the equilibrium 
compositions as function of density and temperature. 

The nonequilibrium flow program is set up in terms of some eighteen sub- 
routines for the various thermodynamic and other functions discussed in Chapters II 
through IV. Its core consists of the standard fourth order variable step Runga- 
Kutta integration routine RKVS, which is described in detail in Ref. 38. Using 
the initial gradients and a simple Euler integration to calculate values for all 
variables at an incremental distance away from x - 0, the computation switches 
over to the Runga-Kutta integration of equations (4.30) through (4.33). Numeri- 
cal problems and running times depend on the case and will be discussed in 
Chapter V . 
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CHAPTER V 

DISCUSSION OF RESULTS 

5.1 Frozen Flow 

The success of the one-strip integral method when applied to the com- 
putation of supersonic flow of a perfect gas ( Y - 1 .4) about circular cones was 
previously demonstrated by South (Ref. 12). He showed that, particularly for 
higher free stream Mach numbers (M^ >3), the results from the integral method 
are in excellent agreement with the charts in Ref. 39. In the present investiga- 
tion, chemically frozen flow was therefore studied only to obtain proper initial 
conditions for the calculation of chemically relaxing flow, as pointed out in 
section 4.3. 

Figures 7 through 12 show the variation of some typical flow field variables 
with cone semi-vertex angle and free stream Mach number, for undissociated air 
in vibrational equilibrium. In Figures 7, 8, 11, and 12, the perfect gas results 
are also given for comparison. In view of previous findings (Ref. 12), the 
difference between the perfect gas results and those of the present investigation 
can be attributed to vibrational equilibrium. It can be observed that for those 
combi nations of free stream velocity and cone angle where the temperature in 
the shock layer becomes high enough to cause significant molecular vibration, 
the shock wave tends to be closer to the body than for the perfect gas case* 
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Simultaneously/ there is a slight increase in the surface velocity, while the 
surface pressure is hardly affected at all . 

5.2 Equilibrium Flow 

Rather than to duplicate the work of Newman (Ref. 36), conical flow 
parameters for air in thermodynamic equilibrium are presented here in order to 
demonstrate the validity of the five component air model for the range of cone 
semi-vertex angles and 'ree stream conditions considered in this work. The use- 
fulness of an equilibrium flow solution in connection with calculations of non- 
equilibrium flow was already discussed in section 4.4. 

Figures 13 through 16 present a comparison of equilibrium flow »esults 
obtained from the one-strip integral method, used in this investigation, with 
results obtained by integrating the Taylor-Maccoll equation (Ref. 21). In 
addition, these figures represent a comparison of the thermodynamic properties 
of the air model used in the present investigation with the thermodynamic tables 
of Blackwell et al ., which were used in Romig's work (Ref. 21). 

Thirteen cases for a variety of cone semi -vertex angles (30 ^ 9 ^45) and 
free stream Mach numbers (10 * s 30) were selected to cover a wide range 

of values of the hypersonic similarity parameter (5 5 M w sin 18) at two 
different free stream pressures. In all cases, the results from the two calculations, 
each using a different air model and a different method of solution, agree quite 


well • 
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Additional eises for conical equilibrium flow were calculated for the 
purpose of comparison with chemically relaxing flow. In particular, it will be 
shown at the end of this chapter, how the state of equilibrium, which is approached 
asymptotically by the nonequilibrium flow, compares with the state of equili- 
brium found when thermodynamic equilibrium is assumed throughout the flow 
field. In the next section, these equilibrium states are, for simplicity, labelled 
“asymptotic equilibrium" and "conical equilibrium", respectively. 

5.3 Nonequilibrium Flow 

Nonequilibrium flow calculations have been carried out for various pur- 
poses, and the free stream conditions were chosen accordingly. Table V gives a 
detailed summary of all cases which were investigated, while Fig. 17 displays, 
for each case, the three independent characteristic parameters (semi -vertex 
angle 9, free stream velocity, and altitude) in a velocity-altitude diagram. 

Cases 1 and 2 were selected principally for the purpose of comparing the 
results from the integral method used in the present investigation, with those of 
the method of characteristics (Ref. 20). Case 3 was computed to check previous 
calculations using a much simpler air model (Ref. 16). Cases 4 through 8 are for 
similar cones at constant altitude but at different free stream velocities. Cases 
9 and 10, together with case 5, were selected to study the effect of different 
altitudes on all flow variables at the same cone angle and flight velocity. The 
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influence of the semi -vertex angle 9 on the flow parameters is studied by com- 
paring cases 10 through 13/ in which both altitude and free stream velocity are 
kept constant. 

Figures 18 and 19 show the species mass fractions on the cone surface for 
cases 1 and 2, respectively. The slight difference in the results near the tip of 
the cone might be due to one or more of the differences in the calculation 
procedures. A small difference in the initial values, for example, which is 
undetectable from the graphs shown in Ref. 20, together with different expressions 
of the equilibrium constants used for the reactions IV and VI (see section 2.5), 
may be one reason for the discrepancy. Another cause may be the fact that, in 
order to use the method of characteristics, a certain small but finite area of 
frozen flow has to be assumed near the tip, while no such assumption is needed 
for the integral method. In any case, the results for the gas composition from the 
two methods are almost identical for x > 4 mm. The temperature and density 
along the cone surface, plotted in figures 20 and 21, show the same excellent 
agreement. 

The variation of the shock wave angle along the length of the cone is 
shown in Fig. 22. The results from the two methods agree insofar as both methods 
predict that the shock wave angle o should decrease. While the characteristics 
solution from Ref. 20 seems to indicate that the shock wave angle cr should 
approach its value for conical equilibrium flow monotomi cal I y from above, the 
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integral method solution indicates a limiting value for cr which is slightly 
higher than that. As a matter of fact, when the integration was carried out to a 
distance of about 70 cm away from the tip for case 2, it was found that the shock 
wave angle has a minimum value (a“ 44.39) between x = 25 mm and x = 75 mm, 
then approaches its asymptotic equilibrium value ( cr = 44.41) from below. From 
x = 23 cm onward none of the variables show any change in the first four sign- 
ificant figures. Apart from this, it is noted that the difference in shock wave 
angle between frozen flow and equilibrium flow is relatively small (only about 
1 .5 degrees in case 1, and even less than that in case 2). 

Cone surface pressures for cases 1 and 2 may be compared in Fig .23. In 
both cases, the difference between the values for frozen flow and those for equili- 
brium flow is less than 2%. The method of characteristics is seen to predict a 
slightly lower surface pressure (lower by about 0.5%) than the integral method. 

Case 1 also was the only case that presented numerical problems. The 
density, the shock wave angle, the pressure and the surface velocity exhibited 
small damped oscillations, of which an example is given in Fig. 23a. It was 
found that a systematic removal of all possible sources for numerical inconsistencies 
in the program reduced the amplitude of these oscillations considerably. A 
further reduction was achieved by finding an optimum size for the initial step 
which was then taken to be A?” 10 for all remaining cases. The fact that 
case 1 has the lowest body surface Mach number, in connection with other 



77 


fir. *ings, indicates that the probable cause for these oscillations lies in a viola- 
tion of the stability criterion for hyperbolic equations. No definite conclusions 
have been reached at the present time since it is felt that this problem requires 
a more detailed analysis. 

The body surface velocity and the surface Mach number for cases 1 and 2 
are shown in figures 24 and 25, respectively. No results from the method of 
characteristics were available for these variables. However, the fact that the 
authors of Ref. 20 present results for wedge flow where u^ was considered to be 
constant, permit the conclusion that they also found that the change in surface 
velocity along the body surface is negligible. 

For both, case 1 and case 2, the flow field has been calculated up to a 
dimensionless distance of 5 = 10.0 from the tip of the cone. This corresponds to 
a physical distance of about 10 mm in case 1, and about 70 mm for case 2- Com- 
putation times were 146 sec for case 1, and 93 sec for case 2. This contrasts with 
computation times of about 1 hour per case on the BRLESC computer (Ref. 20) for 
the method of characteristics calculations. For the integration of case 2 up 
to § - 100.0, the running time was only 252 seconds in this work . 

Figure 26 shows the influence of the free stream velocity on the surface 
species mass fractions for constant altitude and cone semi -vertex angle. Including 
the nitric oxyde mass fractions, which show the lowest values for the highest 
velocity, all variations are as expected. The large difference in the characterise 
tic relaxation lengths is particularly noteworthy. The surface temperatures. 
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shown in Fig* 27, again display the drastic drop near the cone tip, which is 
caused by the strong dissociation gradients in this region. Surface density, shock 
wave anglp, and surface pressure are shown in figures 28 through 30. Their 
limiting values for large ? show only a slight dependence on the free stream 
velocity, in contrast to the gas composition and the surface temperatures. 

Figures 31 through 33, representing a 30° cone at constant free stream 
velocity but at three different altitudes, clearly demonstrate the validity of the 
nonequilibrium scaling law 

p . L s constant 

CD 

which is discussed in detoil in Ref. 40. Here L is some characteristic dimension 
of the flow field or the body. According to the theory, nonequilibrium scaling is 
applicable for that region of the relaxation zone where the reactions are predomi- 
nantly dissociation reactions. As seen from the production rate equations given 
in Chapter II, the dissociation rates are proportional to the density. Once the 
state of equilibrium is approached, the rate of recombination, being proportional 
to the square of the density, becomes of equal importance, and therefore the 
nonequilibrium scaling law should cease to be valid. Figures 31 through 33 show 
that this is precisely so. 

The length scale in these figures is somewhat arbitrary because the reference 
density was chosen for convenience such that the characteristic length 



which should theoretically be the same constant for all cases, was equal to one 
meter for case 10. This means that the x scale is practically identical with the 
§ scale, which permits an easy conversion to the actual physical length in each 
case. 

The reason for the close agreement in all variables shown, in spite of 
differences in the free stream temperature, lies in the fact that, for hypersonic 
conditions, the total enthalpy of the free stream is given mainly in terms of the 
velocity, while the free stream temperature contributes only a negligible amount. 
On the downstream side of the shock wave the situation is almost opposite, and 
therefore the absolute temperature, which controls the reactions, is practically 
the same in the three cases. 

Figures 31 through 33 show that the results, plotted versus x, are identical, 
or at least in very close agreement, up to x = 10. For x > 10, case 10 rapidly 
approaches equilibrium, which is reached near x ” 80. In fact, for x >80, no 
changes in the first four significant figures for the density P and the temperature 
T at the body surface could be observed. The curves for case 9 show a tendency 
towards equilibrium at x = 100, but no such trend can be recognized for case 5- 
It is also noted from figures 32 and 33 that the asymptotic equilibrium values of 
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the surface temperature and ihe surface density do not agree with the respective 
values for conical equilibrium. 

It is concluded that, under otherwise identical conditions, the altitude has 
a decisive influence on the flow field variables. While, for practical purposes, 
the state of equilibrium is reached at a distance of about 50 cm away from the tip 
in case 10 (40 km altitude), this equilibrium state is not yet reached at 150 m 
from the tip in case 5 (80 km altitude). 

The influence of the cone semi-vertex angle 9 on the flow field is shown in 
figures 34 through 37. It is found that, qualitatively, an increase in 9 has the 
same effect as an increase in free stream velocity. Actually, the effect of a 
change in 9 depends very much on the other parameters, namely altitude and 
free stream velocity. It is shown in Fig. 35, for example, that for a variation of 
9 from 20° to 35°, the temperature at the cone tip changes by roughly 8000 °K . 
The values for conical equilibrium indicate that, once equilibrium is reached, 
the temperatures at the body surface still differ by 3000 °K, approximately. 

Figures 38 through 40 display the strong dependency of the characteristic 


relaxation length X on free stream velocity, altitude and semi -vertex angle. 

“3 

For the examples discussed herein, X varies approximately between 10 and 



meters, that is, by five orders of magnitude. 


Finally, Fig. 41 shows a comparison of conical and asymptotic equilibrium 
in a Mollier chart for air in thermodynamic equilibrium. The three equilibrium 
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states shown v ere fixed by their respective values for the enthalpy and the 
density. First, the figure demonstrates that conical equilibrium flow is characte- 
rized by an isentropic compression from the shock wave to the body, since the 
two equilibrium states are located on the same vertical line. Second, it is 
shown that the value of the entropy of the equilibrium state, which is asymptoti- 
cally approached by the nonequilibrium flow, differs from the entropy value for 
conical equilibrium, as discussed in section 4.4. Whether the entropy increase 
shown is mainly due to the difference in the shock wave angle for the two situa- 
tions, or due to the relaxation process cannot be explained unless more detailed 
calculations are carried out. Finally, the close agreement of all variables cal- 
culated for conical equilibrium at the body (e.g. = 5527 °K, ” 1 .255) 

with the corresponding values indicated in the Mollier chart 5500 °K, 

1.253) provides another demonstration for the accuracy of the air model and 
its thermodynamic properties as used in this investigation. 

5.4 Conclusions 

A semi -exact procedure, which uses exact forms of the x-Momentum and 
species continuity equations along the body surface, was shown to yield results 
which are in excellent agreement with those obtained by the method of characte- 
ristics. Although this procedure was also successfully used in the investigation of 
nonequilibrium blunt body flows (Ref. 3,4,5), it was reported to be always 
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unsuccessful for the case of vibrational equilibrium flow past pointed bodies 
(Ref. 12,13). 

The amplitude of some bounded oscillations, which were encountered in one 
case with a relatively low surface Mach number near the cone tip (M^ < 2.5), 
could be reduced to a negligible amount by carefully removing all possible sources 
for inconsistencies in the program. The rapid increase of the surface Mach number 
near the cone tip apparently caused the oscillations to stabilize quickly, and no 
further probl ems were encountered . 

Although it is recognized that the method of characteristics furnishes infor- 
mation about the variables within the shock layer which cannot be obtained by 
the one-strip integral method without additional calculation schemes, it appears 
that running times for the semi -exact procedure used in this investigation are 
drastically shorter than those reported for the method of characteristics. 

The definition of a characteristic relaxation length for dissociation, 

X * (d Z / dx) was found to be extremely useful . Its principal advantage is 

that it transforms all problems to the same length scale, thereby avoiding the need 
to specify integration ranges for each problem individually. The latter process 
could be very time consuming because of the strong variation of X from case to 
case (10 X[m] < 10^, in the examples discussed). 

The application of the nonequilibrium scaling law, . L = constant, to 
cases covering a wide range in altitude clearly demonstrated its validity for the 
type of flow considered in this investigation. 
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It was also demonstrated that, the three parameters, cone semi -vertex 
angle, free stream velocity, and altitude, exert a strong influence on the 
variation of most flow field variables. It is therefore concluded that, for cases 
where the characteristic length is of the same order of magnitude as the length 
of the conical body, large variations of the gas composition, the temperature and 
the density along the streamlines have to be expected in the entire shock layer. 



84 


REFERENCES 

1. Hayes, W. D., Probstein, R. F., Hypersonic Flow Theory , I, Inviscid Flows, 
Academic Press, New York and London, 1966. 

2. Dorodnitsyn, A .A., "A Contribution to the Solution of Mixed Problems of 
Transonic Aerodynamics", Advances in Aeronautical Sciences, Vol . 2, 
Pergamon Press, New York, 1959. 

3. Shih, W . C* L Baron, J . R., Krupp, R. S., and To*le, W. J ., "Nonequil- 
ibrium Blunt Body Flow Using the Method of Integral Relations", MIT Aero- 
physics Laboratory, Technical Report 66, May 1963; also cataloged by ODC 

as AD No. 415 934. 

4. Hermann, R., Thoenes, J., "Hypersonic Flow of Air Past a Circular Cylinder 
With Non -equilibrium Oxygen Dissociation Including Dissociation of the Free 
Stream", paper presented at the Vlth European Aeronautical Congress, Munich, 
Sept. 1-4, 1965; also UARI Research Report No. 28, University of Alabama 
Research Institute, Huntsville, Alabama, September 1965. 

5. Thompson, K. 0-, "Hypersonic Flow of Air Around Blunt Bodies With Finite 
Chemical Reaction Rate", Dissertation, University of Alabama, University, 
Alabama, 1967; also UARI Research Report No. 48, University of Alabama 
Research Institute, Huntsville, Alabama, September 1967. 

6. Allen, H. J., "Hypersonic Aerodynamic Problems of the Future", The High 
Temperature Aspects of Hypersonic Flow , AGARDograph 68, edited by wT C» 
Nelson, The Macmillan Company, New York, 1964. 

7. Stephenson, J. D., "A Technique for Determining Relaxation Times by Free- 
Flight Test of Low-Fi neness-Ratio Cones; With Experimental Results for Air at 
Equilibrium Temperatures up to 3440°K", NASA TN D-327, September 1960. 

8. Sedney, R., South, J. C., and Gerber, N*, "Characteristic Calculation of 
Non -Equilibrium Flows", The High Temperature Aspects of Hypersonic Flow , 
AGARDograph 68, edited by W. C* Nelson, The Macmi I i an Company, New 
York, 1964. 



85 


9. Sedney, R., Gerber, N., “Non -Equilibrium Flow Over a Cone", IAS Paper 
No. 63-71, alsoAIAA J., Vol . 1, No. 11, November 1963. 

10. Capiaux, R., Washington, M., "Non-Equilibrium Flow Pasta Wedge", AIAA 
J., Vol. 1, No. 3, March 1963. 

11. Wood, A. D., Springfield, J. F., and Pal I one, A. J., "Chemical and 
Vibrational Relaxation of an Inviscid Hypersonic Flow", AIAA Paper 63-441, 
alsoAIAA J., Vol. 2, No. 10, October 1964. 

12. South, J. C*, "Application of Dorodni tsyn 's Integral Method to Nonequili- 
brium Flows Over Pointed Bodies", NASA TN D-1942, August 1963. 

13. South, J. C-, "Application of the Method of Integral Relations to Supersonic 
Nonequilibrium Flow Past Wedges and Cones", NASA TR R-205, August 1964. 

14. South, J. C*, Newman, P. A., "Application of *he Method of Integral 
Relations to Real -Gas Flows Past Pointed Bodies", AIAA Paper No. 65-27, 
January 1965. 

15. Newman, P. A., "A Modified Method of Integral Relations for Supersonic 
Nonequilibrium Flow Over a Wedge", NASA TN D-2654, February 1965. 

16. Thoenes, J., "Nonequilibrium Flow Past a Circular Cone Including Freestream 
Dissociation", AIAA J., Vol. 4, No. 8, August 1966; also UARI Research 
Report No. 25, University of Alabama Research Institute, Huntsville, Alabama, 
May 1965. 

m 

17. Thoenes, J ., "Hypersonic Flow of Dissociating Air Past a Circular Cone", 

UARI Research Report No. 26, University of Alabama Research Institute, 
Huntsville, Alabama, January 1966. 

18. Lee, R. $., "A Unified Analysis of Supersonic Nonequilibrium Flow over a 
Wedge: I. Vibrafional Nonequilibrium", IAS Paper 63-40; also AIAA J., 

Vol .2, No. 4, April 1964. 

19. DeJarnette, F. R., "Numerical Solution of Inviscid Hypersonic Flows for 
Nonequilibrium Vibration and Dissociation Using an "Artificial Viscosity" 
Method", Ph.D. Thesis, Virginia Polytechnic Institute, March 1965. 

20. Spurk, J. H., Gerber, N., and Sedney, R., "Characteristic Calculation of 
Flowfields with Chemical Reactions", AIAA J., Vol. 4, No. 1, January 1966. 



86 


21. Romig, Mary F., "Conical Flow Parameters For Air in Dissociation Equilibrium", 
Convair Scientific Research Laboratory, Research Report 7, May 1960. 

22. Herzfeld, K. F., Griffing, V., Hirschfelder, J. O., Curtiss, C. F., Bird, 

R. B., and Spatz, E. L., Fundamental Physics of Gases, Princeton Aero- 
nautical Paperbacks, Vol . 7, Princeton University Press, 1961. 

23. Lee, J. F., Sears, F. W., Turcotte, D. L., Statistical Thermodynamics, 
Addison -Wesley Publishing Co., Inc., Reading, Mass., 1963. 

24. Vincenti, W. G«, Kruger, C. H., Introduction To Physical Gas Dynamics , 

John Wiley and Sons, Inc., New York, 1965. 

25. Wray, K. L., "Chemical Kinetics of High Temperature Air", Hypersonic Flow 
Research , Progress in Astronautics and Rocketry - Vol . 7, edited by F . R. 
Riddell; Academic Press, New York, 1962. 

26. Lin, S. C., Teare, J. D., "A Streamtube Approximation for Calculation of 
Reaction Rates in the Inviscid Flow Field of Hypersonic Objects", Proceedings 
of the Sixth Symposium on Ballistic Missile and Aerospace Technology : 

Reentry , Vol ♦ 4, Academic Press, New York, 1961. 

27. Treanor, C. E-, Marrone, P. V*, "Effect of Dissociation on the Rate of 
Vibrational Relaxation", Cornell Aeronautical Laboratory Report No. QM- 
1626-A-4, Feb. 1962. 

28. Rose, P. H., Teare, D. J., "On Chemical Effects and Radiation in Hyper- 
sonic Aerodynamics", The High Temperature Aspects of Hypersonic Flow, 
AGARDograph 68, edited by W. C. Nelson, The Macmillan Co., New York, 
1964. 

29. Schwartz, M., Green, J., and Rutledge, N. A., ‘ Vector Analysis” , Harper 
and Brothers, New York, 1960. 

30. Wittliff, C. E., Curtis, J. T., "Normal Shock Wave Parameters in Equili- 
brium Air", Cornell Aeronautical Laboratory Report No. CAL- 111, 

November 1961. 

31. Marrone, P. V*, "Normal Shock Waves in Air: Equilibrium Composition and 

Flow Parameters for Velocities from 26,000 to 50,000 ft/sec, " Cornell 
Aeronautical Laboratory Report No. AG - 1729 - A - 2, August 1962. 



32. Sokol nikoff, I. S., Redheffer, R. M*, Mathematics of Physics and Modern 
Engineering, McGraw-Hill Book Company, Inc., New York, 1958. 


87 


33. Busemann, A., "Drucke auf kegelformige Spitzen bei Bewegung mit 
CIberschallgeschwindigkeit", ZAMM, Vol. 9, pp. 496-498, 1929. 

34. Hansen, C. F., "Approximations for the Thermodynamic and Transport 
Properties of High Temperature Air", NASA TR R-50, 1959. 

35. U. S« Standard Atmosphere , 1962, Prepared under sponsorship of NASA, U.S* 
Air Force, U. S. Weather Bureau, U. S. Government Printing Office, 
Washington, D. C., December 1962. 

36. Newman, P. A., "Approximate Calculation of Hypersonic Conical Flow Para- 
meters for Air in Thermodynamic Equilibrium", NASA TN D-2058, January 
1964. 

37. Belotserkovskii, O. M., "Flow With a Detached Shock Wave About a 
Symmetrical Profile", PMM, Vol. 22, No. 2, 1958, pp. 206-219. 

38. UNIVAC 1107 BEEF MATH ROUTINES, Manual No. UP-3984, published by 
the UNIVAC Division of tbe SPERRY RAND Corporation. 

39. Ames Research Staff, "Equations, Tables and Charts for Compressible Flow", 
NACA Report 1135, 1953. 

40. Hall, J. G-, Eschenroeder, A. Q., and Marrone, P. V., "Inviscid Hyper- 
sonic Air-Flows With Coupled Nonequilibrium Processes", IAS Paper, 

No. 62-67, 1962. 


41. Feldman, S., "Hypersonic Gas Dynamic Charts for Equilibrium Air", Avco 
Research Laboratory, Research Report 40, January 1957. 


88 


TABLE I 

GAS MODEL DATA 


Adopted Primary Quantities 

^02*o 

0.21 


( ^N2 ) o 

0.79 


Z “ p /p RT 
0 0 0 0 

1.0 


p o 

1.01325 -JO 5 

CN 

Z 

T 

o 

288.15 

°K 

R* 

8314.32 

J / kmol °K 

N* 

6.02257. 10 26 

kmol ^ 

h 

6.6237 • 10“ 34 

J • sec 

Derived Quantities 

(C 02 ) o 5 ™o 

0.232 918 


(C N2 ) o 

0.767 082 


M 

28.850 335 

kg Amol 

R 

288.187 988 

jAg °K 

P 

0 

1.220 174 

kg/m 3 






z o 


TABLE II 


ATOMIC AND MOLECULAR CONSTANTS 



1) Ref. 35; 2) Ref. 24; 3) Ref. 20; 4) Ref. 34, Values for NO from Ref. 23 








TABLE III 

DISSOCIATION RATE CONSTANTS 



Reaction 


O 2 +M-20+M 


Third Body M 


n 2 , N, NO 


k ^ [ m / kmol • sec ] , Ref . 20 


k d(N2) = '•2 ’0 ,8 T' 1 - 5 «P(- 9 o/T) 

I _ 1 

k d(02) = 3 * C k d(N2) 
k d(0) = h75 * 10 Tk d(N2) 


II N 2 +M z 2N +M 0 2 , O, NO 


III NO+M^N+O+M 0 2 ,N 2 ,0,N,N0 


iv no+o-o 2 +n 



n 2 + O r NO + N 


VI 0 2 +N 2 ;2 NO 


c d (02) 

II 

c d (N2) 

II 

C d(N) 



9.9 -10 17 T“ K5 exp (-9* /T) 


3,03 k d (02) 
15,15 k d (02) 


* 5.2-10 18 T’ L5 exp(-9* o /T) 


- 2.4-10 8 T°* 5 exp(- 19230/ T) 


= 5.0-10 ° exp (-38000/T) 


- 9.1 -10 21 


-2.5 


exp (- 65000 / T) 




















TABLE IV 


EQUILIBRIUM AND RECOMBINATION RATE CONSTANTS 


a 

Equilibrium Constant K 

M 

k = k./K 
r d c 

i 

1.2-10 6 T' 0-5 exp (*0 q 2 A) 
[kmol / m 3 ] 

N 2 ,N,NO 

°2 

O 

.1 - in 12 r 

k r (N2) ° T 

k r(02) =3-0k r(N2) 
k !(C) = '-75-10 9 

ii 

k” - 1.8 -10 4 exp (-^/T) 
[kmol / m 3 ] 

o 2 ,o,no 

N 2 

N 

=55* 10^ 3 T~**^ 
r(02) U 1 

k "(N2) = 3-03 k ”(02) 

k ”(N) = ,5- ' 5kI r I (02) 

hi 

K 1 " - 4.0-10 3 exp(-S^ o /T) 
[kmol / m 3 ] 

o 2 ,n 2 ,o, 

N, NO 

k 111 = 1.3- 10 15 T” 1,5 

r 

IV 

K IV =3.333 * 10‘ 3 T 0,5 exp (-16 1 10A) 
c 


k IV = 7.2 * I0 10 exp (-3 120 A) 

r 

D 

K V *4.5 exp (-37 770 / T) 
c 


k^ = 1.111 -10 10 exp(-230A) 

VI 

K VI - 1 .35 10 3 T' 0 ' 5 exp (-21660A) 
c 


k^ 1 = 6.741 10 18 exp(-43 340A) 


Kj. through from Ref. 25/ kJY through calculated. 



















TABLE V 


TABULATION OF CASES FOR CHEMICALLY RELAXING FLOW 


I 

9*) 

alt.** 

km 

V„*) 

m/sec 

p « . 3 
kg /nr 

T 

°k" 


Comments 

I 

45.43 

31 

gag 

1.287- 10” 2 

273.16 

20 

For comparison with method 

2 

41.07 

II 

■m 

II 

II 

18 

of characteristics (Ref. 20) 

3 

25 

40 

6390.0 

3.996 -10" 3 

250.35 

20 

For comparison with Ref. 16 

4 

30 

80 

11186.0 

1.999- 10 ~ 5 

180.65 

41.4 


5 

ii 

ii 

9350.0 

•i 

II 

34.6 

To stud/ influence of free 

6 

n 

II 

7909.6 

11 

II 

29.3 

stream velocity on flow field 
parameters 

7 

ii 

II 

6750.0 

11 

It 

25 


8 

ii 

ii 

5400.0 

II 

II 

20 


D 

ii 

— 

9350.0 

3.059 -!0" 4 

255.77 

29.1 

To study influence of 

H 

ii 

C S 

ii 

3.996 • I0 " 3 

250.35 

29.4 

altitude (incl . # 5) 

ii 

35 

II 

11 

II 

it 

II 

To study influence of 

12 

25 

II 

H 

II 

II 

11 

semi vertex angle 
(incl . * 10 ) 

13 

20 

II 

n 

II 

It 

II 



*) Parameters which are considered as independent. 
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Fig* 10 Variation of con* surface density with jemf -vertex angle (chem* frozen, 
vlbr. equil., l m * 250 °K). 
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Semi -vertex angle, 0 


Fig. 12 Variation of cone surface velocity with semi -vertex angle (chem. frozen, 
vibr. equil ., T a s 250°K). 



Surface compressibi I ity factor 
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M a sin 0 

Fig* 15 Surface density for equilibrium cone flow (T M = 273.16 K, 
= 1*0 atm)* 
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Integral Method 
Romlg (Ref. 21) 
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Fig. 16 Shock wave angle <7 for equilibrium 

(T_ - 273. 16°K # p =1.0 atm). 
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cone flow 
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Fig* 17 Survey of investigated nonequilibrium flow cases* 



Species moss fractions, 


110 



x [ m m] 


Fig. 18 Cone surface qpecies mass fractions for case 1(8 = 45.43 °, 
V„ = 6638.0 m/sec, J m = 273.16 °K, p m = 0.01 atm). 





Surface temperature 
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Fig. 20 Cone surface temperatures for cases 1 and 2 (T,* = 273. 16 °K 
p w = 0.01 atm). 




Surface density, p, / P 
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Fig. 21 


Cone surface densities for cases 1 and 2 (T^ = 273.16 °K, 
p - 0.01 at»*“l. 

r ao 



Shock wave angle, a Shock wave angle, <7 




b) case 2:8 =41.07°, "5974.0 m/sec. 

Fig. 22 Shock wave angles for cases 1 and 2 O'* = 273.16 °K, = 0.01 atm). 
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b) case 2:0=41.07° , = 5974.0 m/sec. 

Fig. 24 Cone surface velocities for cases 1 and 2 fT^ — 273. 16 °K/ 
p m = 0.01 atm). 





Surface Mach number. 
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Fig • 25 Cone surface Mach numbers for cases 1 and 2 ( “ 273.16 °K, 
P» ~ 0*0^ atm). 



Species mass fractions, 
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5 =x/\ 

Fig* 26 Influence of free stream velocity on cone surface species mass fractions 
( 8 = 30 , altitude 80 km). 





Surface temperature, T . / T 
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Surface density 
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Fig. 28 Influence of free stream velocity on cone surface density (0 = 30 
altitude 80 km). 
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5 =x/x 


Fig. 29 Influence of free stream velocity on shock wave angle (0 = 30°, 
altitude 30 km). 



Fig. 30 Influence of free stream velocity on cone surface pressure 
(0 = 30°, altitude 80 km). 
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? =x/\ 

Fig. 34 Influence of cone semi-vertex angle on surface species mass fractions 
(V„ = 9350.0 m/sec, altitude 40 km). 




Fig. 35 Influence of cone semi -vertex angle on surface temperature 
- 9350.0 m/sec, altitude 40 km). 
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Fig. 36 Influence of cone semi-vertex angle on surface density 
(V ~ 9350.0 m/sec, altitude 40 km). 

CD 
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Fig. 37 Influence of cone semi “vertex angle on surface Mach number 
(V # - 9350*0 m/sec, altitude 40 km). 
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Enthalpy, h / R T 
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Entropy, v'R 


Fig. 41 Comparison of conical and asymptotic epuilibrium for case 10 

(0 = 30°, V — = 9350.0 m/sec, altitude 40 km) in a McMler 2 
diagram (Ref. 41; T_ = 273.16 °K, p - 1.0 atm, p - 1.288 kg/nry 
R = 283. 188 J / kg°K). ° 0 
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APPENDIX A 

SOME EXPRESSIONS RELATED TO THERMODYNAMICS 


Frozen specific heat c^ : 


= R*Jy 02 [ 3 .5 + 


02 


v 2 


+ Y 


N2 


K^r) 


+ Y 


NO 


h <-¥) 


+ 2 5 + V N> } 

Frozen speed of sound/ : 

, c RZT 
J = E 


'f 


Dimensionless functions L, S : 


e 


e v / t 
02 / 1 


(e 9 - /T -,)‘ 


e 


6 n 2 / t 


9 n 2 /t 


s no /t 


4 

-0 


] 

] 


(■ 


0 V /j 2 

NO / 1 




] 


c - RZ 
P 


(A-l) 


(A-2) 


L ~ 


$ = 


qo i 

1 

d x 

(A-3) 

R*T _ 

\ \ s 

, d C. 

1 i 

(A-4) 

V 2 

" m i 

M. d x 

i 
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APPENDIX B 

DIMENSIONLESS FUNCTIONS 0, FOR DIFFERENTIAL SHOCK 
RELATIONS AND INITIAL DERIVATIVES 


2V„ 

COS <7 

V.«in°(V.sina-B)- (2-b) <h $ 

-u" 

m 


T. 

’ (RZ, 

W 

» - 2 c ps ) B + (V.sin o)b e pJ 




V w cos a 
"1 



tan o 






sin o cos a 



where 


b » 1 + 


RZ ro T 
00 00 

(V„sin a ) 2 


(B-2) 


(6-3) 


B 


V 


<V„ sin a f - 2 (h $ - hj 


(B-4) 




^ (sin cr ccs 0 + cos cj sin 



sin a 


sin 

(B~5) 




-1 


^ (sin asin 3- cos a cos 



si n cj cos 3 



The following expressions are valid only at x - 0 for K - 0 : 
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P M V m 

00 00 


(2+j) I 


P u 
s s 


% % 


stn 2 B 


P s 

. 2 „ 

sin P 


+ 2 ( cot 3 + j cot 0 ) ( p w + P s V w n 4 ) (B-7) 

+ j (1 + tan 3 cot 0) (p^ u s + P s ^ ) | 


Cl 


1 


6 7J2 

PcoV^ 


j <P b +p b V" ( Ps +p s u s > 


sin 2 3 


P u v 
$ s s 

—rr 

sin P 


+ j ( 1 + tan 3 cot 0 ) (p^ u s Cl^ + 2 p s V ® u s ° 3 ^ 


»-8) 


+ 2 (cot8 + j cot 9) (p„ u s v $ n 2 + P s v„ v s +P s V. u s n^) 

+ j tan 8 cot 9(P S T„ n, +p„ T s n 2 ) R Z 5 > 


= 


1 


P V 2 

K 0P 00 


( rP U V 
( L si n 2 0 


Ps +P S V S 
sin 2 0 


- j (1 + tan 3 cot 0 ) (p^ v s Cl ^ +P j V # v j Qj + P S V <» u s °4 ) 

- 2 (cot 3 +j cot 0) (p w v 2 Cl 2 +2p $ V,* v s n 4 )l( 2 cot8 + j cot 0) 


> ^ s (P s PVP- T s Q 2 ) 
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APPENDIX C 

FUNCTIONS a.. AND A. IN EQUATIONS (4.30) THROUGH (4.33) 


= U 2 

(C-l) 

11 “b 


12 = " P b c pb 

(C-2) 

21 = RZ b T b + eu b 

(C-3) 

22 =P b[ RZ b * C +S >^3 

(C-4) 

33 = (l + j Scos0)(p„ u s v $ n 2 + P s V„v s n 3 + P s v„ 

« s n 4 ) (C-5) 

44 “ P b “b 

(C-6) 

1 ” y | p s u s [ < 1 + K y s ) tan 3 - ( 1 - K y s ) j 6 S m 0 J 


” p b u b[^ + ^7s^ * an J 


“ P $ v s (l+j 6coi 0) ( 2 + 3 K y^) | 



- % (P« u s Q 2 +p s V ® + J 6cos0 ) — — 

33 


+ 



p, v; 



(C-7) 
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P, Vo 
A = -A 


[ ( ,+e)L b'(x~) %] 

+ —— | <P S + P $ u $ 2 ) [(1+ Ky s ) tan3-(l -Ky s ) j6sin6 J 

' (p b + P b U b 2) [ (,+K>, S )tana+j ‘ sin 9 ^ 

-p u v (1 + j 6cos0) (2 +5 Ky ) 
s s s s 

+ j 6 sin 0 +P s (1 + KyJ ] | 

- e ( 1 + j 6 cos 0 ) F R Z ^ (P $ l o n ] + P„ T s V 


L 

o “1 

+ p u Q« + 2 P V, u n_‘ 
00 s 2 s 00 s 3 


1_L 

j a 33 


A. - 


3 y. 


{ P ,Vs[ ( 


1 + Ky s ) tcm3 - (1-Ky^ j & sin 0 


] 


+ (p, - p ) (2 + j &cos 0 ) - p v ( 1 + j 6 cos 0 ) (2 + 3 Ky ) 

D S S S 5 


+ [ p b~ p s + p b u b +2p S u S 2 0+i 6 COs6) 


K | 


A^ = W.^ + a ( 1 +Ky s ) ( 1 + j 6 cos 0) W Js 


^(c b -S s ) 


P u 

s s 


[(1+Ky s )tan0“(l-Ky s )j 6 sin 0 


- p v (1 + j 6 cos 0)(2 + 3Ky) 
s s s 


+ a (C b - C s ) (1 + j 6 cos 9) (q» u s °2 +P s Va5 


'33 


(08) 


(09) 
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